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ABSTRACT 


Multi-frequency modulation (MFM) has been developed at NPS using both 
quadrature-phase-shift-keved (QPSK) and quadrature-amplitude-modulated 
(QAM) signals with good bit error performance at reasonable signal-to-noise ra- 
tios. Improved performance can be achieved by the introduction of error control 
coding. 

This report documents a Fortran simulation of the implementation of error 
control coding into an MFM communication link with additive white Gaussian 
noise. Four Reed-Solomon codes were incorporated, two for 16-QAM and two 
for 32-QAM modulation schemes. The error control codes used were modified 
from the conventional Reed-Solomon codes in that one information symbol was 
sacrificed to parity in order to use a simplified decoding algorithm which requires 
no iteration and enhances °rror detection capability. Bit error rates as a function 
of SNR and £,/N) were analyzed, and bit error performance was Weighed against 


reduction in information rate to determine the value of the codes. 
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THESIS DISCLAIMER 


The reader is cautioned that computer programs developed in this research 
may not have been exercised for all cases of interest. While every effort has been 
made, Within the time available, to ensure that the programs are free of compu- 
tational and logic errors, they cannot be considered validated. Any application 


of these programs Without additional verification is at the risk of the user. 
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I. INTRODUCTION 


Many modern communication methods which merge communications and 
computer technology require efficient and reliable transfer of digital data from 
an information source to a destination. This requirement is not always an casy 
one to achieve. Because this transmission is through a physical medium which is 
subject to various tvpes of noise, distortion and interference. the data output from 
the channel mav differ from its input. The demand for reliable communication 
necessitates the detection and possibly the correction of those errors resulting 
from impaired transmission. These errors can be detected and corrected using the 
principles of error control coding. 

Error control coding can be implemented in a communication svstem usually 
at the cost of data transmission rate. Generally, the code that carrects more er- 
rors Wil require more parity symbols. reducing the data rate. On the other hand. 
a code Which corrects more errors Will be more reliable in a high noise environ- 
ment. The balance between data transmissien rate and coding gain poses an in- 
teresting problem. one which will be addressed in this thesis. 

The subject of this thesis is the simulation and analvsis of various error con- 
trol codes in a bandpass Multi-Frequency Modulation (MFM) communication 
link with Additive White Gaussian Noise (AWGN). MFM is a signal modulation 
format Which is readily adaptable to many digital communication link scenarivs. 
Its efficient use of the frequency spectrum allows very high data transmission 
rates. But like all communication systems. MFM is subject to errors. 

There are four error control codes used in the simulation, all of which are 
Reed-Svlomon. Two error correction detection capability codes (single error 
correcting double error detecting (SEC DED) and double error correcting triple 
error detecting (DEC TED)) in combination with three different modulation 
schemes (one 16-QAM and two 32-QAM constellations) provide performance 


data for the four codes. The effectiveness of implemented codes, in terms of bit 





error rates. coding gains obtained versus the uncoded link, and code rates are the 
areas of study for this report. 

Chapter Two is a brief overview of MFM as it pertains to a 16-20 kHz 
bandpass channel. the communication link under simulation. Chapter Three is 
a discussion of Reed-Solomon error control coding in general, including the en- 
coding and decoding processes and other parameters necessary to define that 
class of codes. Also included is a descripuon of how the codes implemented difter 
from conventional Reed-Solomon codes. Chapter Four describes the link simu- 
lation program and supporting programs and files, channel coding and decoding, 
how error control was implemented into the MFM environment, and the results 
obtained using the various codes compared to the uncoded channels. Chapter 
Five contains the conclusions and recom aendations of possible wavs to achieve 


Improved performance. 





Il. THEORY OF MULTI-FREQUENCY MODULATION 


A. BACKGROUND 

Most modulation methods for bandpass channels hold the digital signal in- 
formation in the arr litude and or phase of a carrier (e.g., Phase Shift Keving 
(PSK), Quadrature Amplitud2 Modulation (QAM)). Information may also be 
transmitted by the presence of one of two carrier frequencies (Frequency Shift 
Keving (FSK)). MFM is a variation on these methods. The signals are directiy 
encoded. modulated. demodulated. and decoded using digital signal processing 
techniques contained in the software of the host computer. 

In MFM. signals are grouped into “packets” which are arbitrarily located in 
the frequency spect:um, and consists of one or more bauds. An MFM packet is 
shown in Figure 1, on a disolay of frequency and time. This figure shows the 
nature of MFM: a structured block of time and frequency slots. 

This structure makes MFM a technique that is ideail* suited to ailow the host 
computer to absorb the functions of modulation and multiplexing, demultiplexing 
and demodulation {Ref. 1: p. 3]. The packet consists of " bauds, each holding 
information on K tones. The baud length is AT during which K discrete tones 
are transmitted. In MFPSk, the phase of each tone contains the information, 
and in MFFQAM, both the phase and magnitude of cich tone contain the infor- 
mation. These Lk subsignals form an orthogonal signal set. Each of these sub- 
signals may be independently modulated with both phase and amplitude 


information. 


B. GENERATION AND DEMODULATION 

Generation and demodulation of MFM is based on the symmetry properties 
of the iast Fouricr Transform (FFT). That is, the FFT of a real-time signal 
corresponds to a conjugate symmetric frequency spectrum, and the inverse FFT 
of a conjugate svmmcetric frequency spectrum corresponds to a real-time signal. 
This means that. if oniyv half of the frequency spectrum is used for in-phase and 


quadrature information, and the other half is loaded with its complex conjugate 
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Figure 1. NIFM Signal Packet: (after Ref. 1: p.3) 


image, then an inverse FFT of this spectrum will result in-a real-time signal. 
Conversely, the FFT of this real-time signal will result in the same frequency 
spectrum. Thus, the MFM transmitter is a device which takes an inverse FFT, 
and the MFM receiver is a device which takes a FFT (with appropriate D A and 
A D conversion). 

In a bandpass channel, the usable portion of the frequency spectrum for in- 
formation is further reduced. Current MFM applications under study at Naval 
Postgraduafe School deal with a 16-20 kHz bandpass channel, with a D.A and 
A.D sampling frequency of 61440 Hz. The corresponding parameters for the 
MFM packets are shown in Table 1. 

In the communication link under study, the baud type(s) transmitted in any 
give~ ,.-ket is known by the receiver. Because the tones containing information 


occur in a known band of the frequency spectrum, only those harmonic numbers 








of the FFT corresponding w tnis band need be analyzed. This function is similar 
to bandpass filtering. but noise within the pass band is still present, leading to 
potential errors. 


Table 1. MFM PACKET PARAMETERS (BANDPASS 16-20 KHZ) 


BAUD TYPE 4 
Space Between 120 

Tones Af (Hz) 
Lower Harmonic A 7 


ee ee ee ee 


Upper Harmonic * 


Upper Tone (Hz) | 19920 | 19920 19920 
Number of Tones 156 
With Information = 

FFT Size 56 512 4 5 























II. REED-SOLOMON CODES 


A. INTRODUCTION 
A binary block error control code with block length n which contains 2* code 
words is called a linear (n,k) code if and only if those 2‘ code words form a k-di- 
mensional subspace of the vector space of all the n-tuples over the Galois Field 
GF(2). A cyclic code is a linear block code in which every cyclic shift of a code 
word (left or right) is also a code word. Cyclic codes are useful because the en- 
coding and syndrome calculation processes can be accomplished easily through 
shift registers with feedback connections and because their inherent algebraic 
structure makes possible the use of several practical methods for decoding 
{Ref. 2: p. 85]. These codes are normally expressed as polynomials. A large class 
of powerful cyclic codes are those known as BCH codes. Their importance in 
coding is due to the fact that they are capable of correcting all random patterns 
of ¢ errors with a simple decoding algorithm that is easily implemented [Ref. 3: 
p. 167}. . 
Reed-Solomon codes form a special subclass of non-binary BCH codes. A 
conventional (7.k) t-error correcting Reed-Solomon code which is comprised of 
elements of the Galois Field GF(q) has the following structure: 
® code word block length n=q-—1 
® number af parity symbols n-k=2t 


® minimum distance Aan = 2t4+ 1 


For this design, g is taken to be a power of two (i.e., g = 2” ). Appendix A dis- 


plays a binary representation of the elements of GF(2*) and GF(2°). 


B. ENCODING 
If « is a primitive element of GF(2”), then the generator polynomial of this 


(n,k) Reed-Solomon code is given by 











G(x) = (x + a(x + 0”)... (x + 074) 
21-1 20 (1) 


2 
= Po Pe kyx FH Box" He Bon * + x 
The generator polynomial roots can also be shifted to save circuit components 


[Ref. 4: p. 174). In this case, G(x) takes the form 
G(x) = (x + (x + TY) (x + 0 FA), (2) 


where jp is any integer. The code generated by G(x) consists of all polynomials 
with coefficients from GF(2”), of degree n — 1 or less, and which are multiples 
of G(x). The message to be encoded is a polynomial, U(x), of degree 
k ~1=n—2t-—1 or less. In non-systematic form, the defreebonaine codeword 


polynomial is 
V(x) = G(x)U(x) = vg + ya + vax? Ss ion (3) 


This is clearly a multiple of G(x) . 

Encoding in systematic form requires more computation. To structure the 
codewords with the information svmbols preceeded by the parity symbols, U(x) 
must first be pre-multiplied by x2. Then x#U(x) is divided by G(x) , which yields 


a quotient Q(x) and a remainder B(x) of degree 2¢ — | or less: 
x(x) = Q(X)G(x) + BOX) (4) 


If B(x) is added to x?"U(x), then the result is a multiple of G(x) which is the sys- 


tematic codeword corresponding to U(x): 





x"U(x) + B(x) = Q(x)G(x) (5) 


This polynomial division can be implemented in a shift register, feedback circuit 


as shown in Figure 2. 


Message ¥7~* u(X) Q Code word 


Parity-check 
digits 





Figure 2. Systematic Encoding Circuit for an (n, k) Cyclic Code: [Ref. 2: p. 95] 


Using the circuit in Figure 2, an (n,k) linear cyclic code with generator 


polynomial G(x) as defined in equation (1) or equation (2) is encoded as follows: 


1. With the gate turned on, the k information symbols uw, u, ..., 4; are shifted 
into the communication channel and at the same time shifted into the the 
encoding circuit. This automatically has the effect of pre-multiplying U(x) 
by x". 


2. After k clock cycles, the entire message has been shifted into the channel 
and the n—k shift registers contain the elements of the remainder 
polynomial. 


3. The gate is opened to break the feedback connections and the switch is 
shifted to transfer parity symbols. 





4. The parity symbols are clocked out onto the channel. These »—k parity 
symbols together with the k information symbols form the a-length 
codeword. 


Reed-Solomon codes in systematic form use the same encoding circuit. 


C. DECODING 

Decoding Reed-Solomon codes can be performed using a table-lookup 
method or using iterative algorithms (Berlekamp, Forney, Euclid, etc.) (Refs. 4, 
5}. For binary BCH codes, these algorithms determine the location of errors 
(using the error location polynomial). For non-binary codes, as in the case of 
Reed-Solomon codes, once the error location polynomial is found, the magnitude 
of the error(s) must also be found in order to correct the received signal block. 

Because every codeword }(x) is divisible by the generator polynomial G(x) 
and the roots of G(x) are «, «*, ..., «2, these are also the roots of V(x). As such, 


the following equation holds for 1 <i< 2t: 
on Pan Se adel 2 (n—1)i _ 
P(x)=votvya tye” t+ yo =0 (6) 


or in vector form, 











By forming ann by n—A = 2t matrix whose columns are the vectors of powers 


of g, §=1.2....2f in the previous equation, any codeword vector will be 
orthogonal to that matrix. The matrix thus formed is the transpose of the parity 
check matrix for a conventional (7.4) Reed-Solomon code. In other words, the 


parity check matrix for an (.k) Reed-Solomon code is given by: 
2 ul 
l uae vA 


fy) “Wy a ey" 
1 (z") (x3)? cee yr 


cas ae 


The syndrome of a received polynomial is found from multiplying the re- 


ceived polynomial (in vector form) by the transpose of the parity check matrix: 





s=reH'=(vy+e)+H' =eeH!, (9) 


where s is the syndrome vector, r is the received vector, v is the transmitted 
codeword vector, e is the error vector, and H is the parity check matrix. (Recall 
that any codeword vector v is orthogonal to H'.) Syndrome component s, can 
also be calculated by evaluating the received polynomial, A(x), at «', for 
i= 1,2.....22. Once the syndrome components are found, they are used in the it- 


erative calculations needed to find the error location polynomial, ¢{x), where 


c(x) 


(1 + Byxya + Pox) oe (i + Bx) 


(10) 


< , 
Gg toxttow te +0.x° 


Here. the roots of a(x) (fyi. Bs', ... Bo!) are the error location numbers, and v is 


the total number of errors in the received codeword (for v <¢). Then defining 


L(x) = 1 + (5) + O))X + (52 + 6)5; + 2)X + 


| y (11) 
+ (s, + G50) + 25,_9 Fo ta,)N. 
the error magnitude at location f, = # is [Ref. 2: p. 175]: 

Z(B7') : 
ae a, (12) 

-1 

(1 + Bf; ) 

i=] uAl 


This demonstrate. the ease of implementation of the iterative algorithms: the 
3 


disadvantage is that the iterative solution does not provide high-speed decoding. 





On the other hand. a table-lookup method can provide high-speed decoding, but 
requires too much memory or overly complicated logic circuitry even for moder- 
ate code lengths [Ref. 5° p. 1]. However, some modifications to Reed-Solomon 


codes can allow a faster decoding procedure. 


D. MODIFIED REED-SOLOMON CODES 
There are several ways to modify Reed-Solomon codes to fit specific applica- 
tions. The block length can be extended or shortened. and the number of parity 
or information symbols can be varied. The (2"—1, 2" —2r — 2) codes utilized in 
the simulation carry fewer information symbols (more parity) than the conven- 
tional codes. reducing the code rate. The benefits lie in the enhanced error 
correction detection capability and speed of decoding. 
I. The (2"— 1,2” — 4) Reed-Solomon Codes 
By ieplacing one of the information symbols in the conventional 
(2"—1, 2" —3) code with a parity symbol. the single-error correcting code can 
also detect all double errors. This modified SEC DED code has the following 
characteristics: 
e codeword block length n=q-—1 
@ number of parity svmbols n-k=3 
® minimum distance din = 4 
Generation of this code is similar to conventional codes, except that now 


the generator polynomial has an additional factor: 


(x + I)(x + a)(x + x”) 
{" tally 4 glx? 4x3, for the (15,12) code (13) 
ats 


G(x) 


3 
wd 


D) » 
wx tal x? + x?) for the (31.28) code 











The encoding circuit for these codes is still that of Figure 2. The major difference 
between these modified codes and the conventional codes lies in the decoding. 

Reference 6 details a fast decoding technique for an extended SEC.DED 
Reed-Solomon code based solely on the syndrome components (i.e., iterative al- 
gorithms are not required). The same decoding technique applies to the modified 
codes used in the simulation, since the modified codes are shortened versions of 
the extended codes. For the simulation, the blocklength of the extended codes is 
shortened by three symbols to better match the number of tones per baud type, 
and to facilitate encoding. 

The extended codes of Reference 6 are formed by the addition of a three 
by three identitv matrix to the parity check matrix of the code generated by the 
generator polynomial of equation (13). This parity check matrix is then short- 
ened by deleting the three right-most columns for use in the (27-1, 2” —4) sim- 


ulation codes: 


1oo011 4] 
H=|0 101 2% (x)? (x)"~4 (14) 
C.0 DL yk ey 


The decoding algorithm for these codes is as follows: 


1. Calculate the syndrome components (in vector form, $ = (5, 5;, 5;)). If all are 
zero, decide no errors occurred, 


2. If a single error occurs in one of the first three positions of the received vec- 
tor (Which cerresponds to the identity matrix positions of H). then the syn- 
drome has onlv one non-zero element. The error has magnitude equal to the 
non-zero svndrome component and location corresponding to the svndrome 
component number. 


3. A single error in any other position satisfies 





ll 

i 

R 
_—~" 
— 
nN 
~~ 


With the error in position ¢ of magnitude 5. 
4. If a single syndrome component is zero, or if 


S| S45 
= # a (16) 


then decide that at least two symbol errors occurred. 


5. Once errors are corrected. recalculate the syndrome to ensure all components 
are zero. If not, decide that two or more symbol errors occurred which 
yielded the same syndrome as a single error. (Note: this step was not in- 
cluded in Ref. 6 but was found necessary to detect all double errors.) 


AS an example. suppose that for the (15.12) code a single error of magnitude «? 
occurred in the fifth position of the received block r. Then the error vector is given 
by e = (0.0.0.0. 2’.0.....0). and the svndrome is s=reH! =eeH! = (2°, x!!, g!3), 
Then s, & = s,s, = 4° = 4°}, corresponding to an error of magnitude 5 = «° in the 
fifth position, 
2. The (2”— 1,2" — 6) Reed-Solomon Codes 

The (2”— 1,2” — 6) codes are obtained by making the same alteration to 
a conventional double-error correcting code as was the case for the SEC. DED 
code: sacrifice one information symbol to parity. This now allows the detection 
of all triple errors. and as will be shown later, the correction of some triple burst 


errors. The characteristics of the DEC:TED codes are: 
¢ code word block length n=q-1 
® number of parity symbols n-k=§ 


® minimum distance dan = 6 


The generator polynomial for this code also has an additional factor over 


the conventional code. and using shifted roots is given by: 





G(x) = (x + a? V(x + a (x + I(x + a(x + a”) 


{) tatx tot? 4 ott txt 4 x9) for the (15,10) code (17) 


1D bh tl x? te ol xP + abet 4 x, for the (31,26) code 


Once again, these codes are encoded in systematic form using the circuit of Figure 
2. Decoding is, however, somewhat more involved than for the SEC’ DED codes. 

Reference 5 gives a detailed description and proof of a fast-decoding 
technique for these codes which is again solely based on the syndrome of the re- 
ceived signal, The parity check matrix of these codes has an additional row over 


the conventional codes and is given by: 


H=|! 1 ! I (18) 
ie er (x) 
I x (3 e (x7)"7! 


The steps of the decoding algorithm are as follows: 


1. Compute the svndrome (in vector form, $ = (5_., 5_4s Sp. 5; 5:)). If all compo- 
nents are zero, decide that no errors occurred. 


ho 


If more than two components of the syndrome are Zero, decide that at least 
three symbols are in error. 


3. Compute y,, 72, and y;, where 


Vy = 51S2 + $15 (19) 
ee 2 5 
2= S25_9 + So (20) 

¥3 = SpS, + $25_}- (21) 








iY 
? I] : ‘ 
If +, = 72 = 7; = 0 then calculate 7’ =——. Correct a single error at location 
i of magnitude 5. . 
4. If s).y.,and +, are not all zero but at least one of them equals zero, decide 
that at least three symbols are in error. 


, 


5. If none of the +,’ are zero, compute k =c.6? and the trace of k. where 
Faeyr «ce = sy 4 and the trace of k, ancelement of GF(2"), 1s defined ‘as 


m—] 
Trik)= ) k?, (22) 
i=0 
If Tr(k) = 1 , decide that at least three svmbol errors occurred. 


6. If Tr(A) = 0. assume that a double swmbol error occurred and sulve 
yet by +ec=0 to find the roots z and #. Then the first error is in the 
i position with magnitude 

—— (sy% + 5) 

(asses b 


and the second error is in the y* position with magnitude 


C2 = Sg t+ Cy. (24) 


~i 


. Onee errors are corrected. recaiculate the svndrome to ensure all components 
are zero. Tf not, decide that at least three symbol errors occurred which 
Vielded the same svndrome as a double error. (Note: this step was not in- 
cluded in Ref. 3S but was found necessary to detect all triple errors.) 


Reference 5 does not examine the code’s capability of triple error cor- 
rection, but the correction of certain triple burst-errors is possible. Suppose that 
three errors of magnitude e, occur in consecutive positions of the received signal, 
the frst at posiuion 7, §= 0.1... — 1. (Note that wrap-around errors can also 
be corrected, e.g.. errors in the (#—1)",O*% and 1" positions. ) The error 


polynomial £(x) is of the form 


— . ; 
E(x) = ex) + ext! + ex'*?, (25) 





or in vector form, e = (0.0,...,0.¢,.e,, @,.0,....0). Multiplying the error vector by the 


transpose of .he parity check matrix gives a syndrome vector of the form 


$_y c(t + git) ae 2420 
Sa ate ee Oy a Oh 
ius i (26) 
1 ey(a + ait! + at?) 
$2 e, (a! 4 FD) o2li#2)) 


From these relations, it is clear that if this type of triple error occurs, then the 


error locations and magnitude can be calculated solely from the syndrome if the 


following simplified conditions exist: 


S_2 ea (1+ a? +074) 

S_y ea ban + a*) 

So] = ej (27) 
5 ea'(l +a +07) 

53 e,a"(1 + a? + a’) 


Note that the addition of elements of GF(2”) vary with the value of mm, so the 
exact relations between syndrome components and error locations vary with the 


code used. As an example. for the (15,10) code with elements from GF(2‘), the 


relations will be 


Sy Bee) eat? 
S_y ex (a8) eae! 
s=}s {=| ¢ [=] e (28) 
7 ¢,2'(a!) elt 
$2 ey2""(a*) eat 


Then to check for this type of error in the (15,10) code, the syndrome components 


must satisfv s_s; = (s,)°z* and s_.s, = (s)*2% The error magnitude is then the 
Value of sy and the starting location can be obtained from any of the other syn- 
drome components, say by dividing s, by s)z!° and taking the exponent of «. 
This capability allows correction of an additional 240 errors [15 triplets x 16 ele- 


ments in GF(2+*)] for the (15,10) code and 992 errors [31 triplets x 32 elements in 


GF(2°*)} for the (31.26) code. 








IV. SIMULATION AND RESULTS 


A. BACKGROUND 

The MFM link is simulated in the Fortran-77 program “MFMLINK”. For a 
given baudtype and code (uncoded 16-QAM, uncoded 32-QAM, and (15,10), 
(15,12), (31,26) and (31,28) Reed-Solomon codes), the program chooses random 
data from GF(2'*) for 16-QAM or from GF(2°) for 32-QAM. Uncoded QPSK was 
also included for comparison with the QAM results. 

For the cases with error control, this data is encoded into Reed-Solomon code 
word blocks. For each baud, the elements of GF(2”) are ‘en transformed into 
the appropriate QAM symbol with each tone representing one element of 
GF(2”), loaded into the frequency spectrum (with conjugate symmetry) and 
transmitted with the inverse-FFT. At this point (in the time domain), white 
Gaussian noise is added. The FFT is then taken, and each FFT sample within 
the information band (from tones k, to k,) is transformed into an element of 
GF(2”). 

For the coded cases, these k,—k,+I1 symbols are divided into code word 
blocks, and each block is decoded separately. The output of the link is then 
compared to the input for error calculation. A block diagram of the link is shown 
in Figure 3. A summary of the Reed-Solomon codes used is shown in Table 2, 
with code rate R=k/n and spectral efficiency 7 = Rye (Nun iS the uncoded 


MF-QAM spectral efficiency). 








Additive 
White 
Gaussian 
Noise 


MFM Transmitter MFM Receiver 


Reed-Solomon Load Spectrum/ FFT/ Reed-Solomon Estimated 
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Figure 3. Simulation MFM Link 


Table 2. SUMMARY OF CODES IN SIMULATION 
Block eo (22) 
suai Symbols 


Menem eae 
SE AEROS TEC NOS NED RE 


: Enor 
SEC/DED | DEC/TED* SEC/DED DEC/TED* 

















* can correct some triple burst errors 


B. CHANNEL CODING AND DECODING 
The 16-QAM constellation consists of two magnitude rings and eight phases. 
The magnitude rings are of amplitudes 1V and 2, with the value of Y variable. 


Data from simulation runs was collected with Y = 100. The phases employed 


were offset +22.5° from the axes. Figure 4 is a diagram of the 16-QAM 








constellation and Appendix A demonstrates the one-to-one relation between 
Galois field elements and QAM constellation elements for both 16-QAM and 
32-QAM. Decoding for this constellation is accomplished for each sample in the 


information band of the FFT as follows: 


1. Determine the quadrant of the received complex signal by the sign of its real 
and imaginary parts. 


2. Determine which of the two phases in the quadrant is closest to the signal. 
This is done by comparing the absolute value of the real pait to that of the 
imaginary part. If that of the real part is greater, choose the phase closest to 
the real axis. If not, choose the phase closest io the imaginary axis. 


3. Finally, estimate the correct magnitude: compare the received magnitude 
with 1.5X. If the received magnitude is greater than 1.5X, assume that the 
received signal is actualiy the constellation element on the outer ring. If not, 
assume it is the clement on the inner ring. 


Two 32-QAM constellations were used in the simulation. The first is Gray 
coded and consists of four magnitude rings and eight phases. The magnitude 
rings are of amplitudes 1X, 2X,3X and 4X, once again with XY variable. Data 
for 32-QAM was collected with X = 20. The eight phases are the same as for 
16-QAM, + 22.5° offset from the axes. Figure 5 is a diagram of this Gray-coded 
constellation. The second 32-QAM constellation is identical to the first except 
that the rings of magnitude 2Y and 4X are shifted clockwise 22.5’ to increase the 
minimum distance between constellation elements. Figure 6 is a diagram of the 
modified 32-QAM _ constellation. Channel decoding for both 32-QAM 
constellations uses a different algorithm than 16-QAM, one which finds that ele- 
ment of the constellation that is closest in distance to the received signal (a form 


of maximum-likelihood decoding). It is described as follows, with r, being a re- 


ceived complex sample in the information band of the FFT and q,, i=0,],...,31 


being the elements of the 32-QAM constellation: 








Figure 4. 16-QAM Constellation 


to 


. Let Yuin = J (Re(r,) — Re(qo))? + (Im(r,) — Im(q))?_, the length of a line seg- 


Vy Jo 
ment drawn between r, and q) on the real-imaginary axis. 


- COMPAL Yin tO Yj = V (Re(r,) — Re(q,))? + (Im(r,) — Im(q,))? for each suc- 


cessive element of the constellation. If y, 1s less than }'nia» replace the value 
Of Yan With the value of y,, store the identity of the constellation element 
Which gave the smaller value (denoted by the subscript 7), and try the next 
element. 


. After comparison to all elements of the constellation, ¢ identifies the closest 


element. Assume that the element / is the actual symbol. 


The QPSK constellation consists of four phases: one each at +45° and 


+ 135°. Data was collected with these points having magnitude ten. Channel de- 


coding consists of determining the quadrant of the received samples based on the 





Figure 5. Gray Coded 32-QAM Constellation 


signs of their real and imaginary parts. The received phase is assumed to be that 


which lies within the reccived quadrant. 


C. SUPPORTING FILES AND PROGRAMS 

Prior to Writing an MFM link program, two data files were created ("“GF16 
DATA” and “GF32 DATA”) which contain the binary representations of the el- 
ements of GF(2*) and GF(2°), respectively. These files form the basis of the 
Reed-Solomon computations in the link. In addition, three data files ("QAM16”", 
“QAM32” and “QAM32A’) were created which outlined the QAM constellations 


used, containing magnitude and phase information on cach element. 





Figure 6. Modified 32-QAM Constellation 


Several other programs were written to support the operation of MFMLINK. 
The program “GALOIS” reads the si-bit representations of a Galois field and 
outputs the GF(2”) addition table. GALOIS adds each element bit-by-bit 
(modulo-2) and then compares the result against each element of the input data 
file to determine which element of GF(2”) it corresponds to. This addition table 
is Output to the data file “"GFAI6” or “GFA32” as. appropriate. The program 
“BITDIFF” reads in the m-bit representations of GF(2”) and compares each ele- 
ment bit-by-bit to determine the number of places in which cach differs for use 


in calculating the number of bit errors. The output of BITDIFF is a GF(2”) 





“difference” table, either “GFD16” or “GFD32”. The 16-element difference table 
for the QPSK symbols was entered separately in the file “QPSKDI”. “TRACE” 
reads in the Galois field addition tables created by “GALOIS” and computes the 
trace of each element (see eqn. 20). Finally, “"QAMCOM” reads in the magnitude 
and phase information of each QAM constellation and calculates the corre- 
sponding real and imaginary components. This output is placed in data file 
"QAMRII6"”, “QAMRI32 or “QAMRI32A”", for use in channel encoding and 
decoding. The real and imaginary components for the QPSK constellation were 
entered separately in the data file “QPSKRI”. 

The outputs of MFMLINK are written into the data files “RSRSLTS” and 
another file with a variable name, depending on the coding scheme used. 
RSRSLTS contains information on code chosen, total symbol and bit errors, 
number of bauds run, and (if desired) the actual random data symbols input, 
transmitted code words, received signals, and corrected code words. The file with 
the variable name contains only the SNRs used in the run and resulting bit error 
rates for use in plotting, and is named according to the code and baud type of the 
run (e.g., “COBT1” corresponds to code “0” (uncoded 16-QAM) and baud type 
“E"): 

D. MFMLINK 

MFMLINK is subroutine based, for euse in altering parameters or algo- 
rithms. One of the difficulties in adapting a finite Galois field into a software 
simulation is the use of non-binary field arithmetic, the addition and multipli- 
cation of the powers of primitive element «. In the program (and supporting files), 


this was handled by identifying ficld elements in terms of its power of « plus one, 
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except for the zero element which is always “0” (i.e.,0-0, l=a° > 1, a! 42, 
etc.). With this method of element identification for GF(2”), field 
multiplication division of elements other than zero is carried out by 
addition, subtraction modulo n = 2” —~ 1. 

A program run starts with the user choosing the type of coding scheme and 
baud, and, if error contro} is utilized, whether code words with more errors de- 
tected than the code can correct should be erased. SNRs and the number of 
bauds for each SNR are fixed in the program and must be edited in order to 
change. Once the user inputs these parameters, subroutine "BAUDS"” reads the 
FFT size and its logarithm base two, k,, k,, and the number of code words per 
baud (when error control is used). Subroutine “DATA” then reads in the Galois 
addition and difference tables, element traces, and QAM constellation charac- 
teristics. Once this background data is read, the run commences. 

If an uncoded case is chosen, subroutine “NOCODE” provides random 
Galois field elements as data and the corresponding QAM or QPSK constellation 
clements, and loads the inverse-FFT array for transmission. For coded cases, the 
main program determines the random Galois field data for a code word, which 
is then encoded in Reed-Solomon systematic form (with subroutine “ENCODE’”). 
This process is repeated until the baud is filled with Reed-Solomon code words, 
with each tone representing one element of the Galois field in use. The main 
program then loads the corresponding inverse-F FT array for transmission. 

Subroutine “BITREV” places the array in bit-reversed order for use in the 


subroutine “INVFFT” which takes the inverse-FFT of the transmission array 





using the decimation-in-time algorithm. Both BITREV and INVFFT were copied 
from the digital signal processing files of Professor Robert D. Strum at NPS. 
Once a baud has been transmitted, noise is added to each sample of the time 
domain signal. First, signal power is anleulatea from the sum of the squares of the 
time domain signal samples (a real-valued sequence). Then noise power (or var- 


lance, Since the noise is fixed to have zero mean) is calculated from 


Spine nee Be 
Ba 1) 1 9(S.¥R-n)/10.0 (29) 
where o2 is the noise power, N, is the FFT size, 7, is the inverse-FFT array 


(transmitted signal), and 7 is given by 


1 = 10 logo (30) 


2 
ky —ky 


The quantity 7 corrects the value of SNR to be the ratio of signal power to noise 
power in the frequency band from k, to k,, rather than having the signal power 
concentrated in this frequency interval and noise power spread uniformly 
throughout the spectrum. Once the noise variance is calculated, subroutine 
“GAUSS” calculates random samples of AWGN with zero mean and variance 
o2 to add to each sample of the time domain signal, forming the received signal 
samples. 


The subroutines BITREV and “FFT” (FFT is also copied from Professor 


Strum’'s files) function as the MFM receiver. Each reccived complex sample in the 








information band from k, to k, is then channel decoded into an element of the 


appropriate Galois field by the subroutines “WHAT16” (for 16-QAM), 
“WHAT32” (for 32-QAM) or “WHATQP” (for QPSK). For uncoded cases, these 
decoded symbols are compared with the input symbols to determine symbol and 
bit errors. For runs with error control, after channel decoding each estimated 
code word is then decoded using the appropriate algorithm of Chapter III. First, 
each estimated codeword’s syndrome is calculated (with subroutine “SYNDRO’). 
If any errors are detected, either subroutine “DECONE” (for SEC codes) or 
“DECTWO?” (for DEC codes) are called to correct errors, if possible. If code word 
erasure was called for by the user, these decoding subroutines call “RUBOUT” 
to erase code words that cannot be corrected. After Reed-Solomon decoding, the 
input code words are compared to the corrected output to determine the number 
of symbol and bit errors. Finally, subroutine “OUTPUT” sends pertinent output 


information to the data file RSRSLTS. 


E. RESULTS 

The primary results of the simulation are displayed in the plots of bit error 
rate versus SNR and bit error rate versus E;/N, . At lower values of SNR (or 
E,/No), accurate results were obtained by performing runs with 10,000-20,000 
bits. At higher values of SNR, as the number of errors diminished, over 100,000 
bits were required. The results of the simulation runs are plotted in the eight 
graphs included in the following pages. The graphs are in pairs and reflect the 
bit error rates as a function of either SNR or E;/Ny . The plots with £,/N, are 


derived from those with SNR based on the following relationship: 
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m, (31) 


where SNR, is the wide band SNR (that used as an input to the program) and 
mis the number of bits per tone for a particular modulation scheme. The calcu- 
lation of E,/N» includes the energy in all bits carried on a tone, both information 
and parity, since the bit error calculations include errors in both data and parity 
symbols. Figure 7 is a graph of the uncoded modulation schemes’ performances 
on the link in terms of Bit Error Rate vs. SNR_~ and also shows the same per- 
formance in terms of £,/\, . Figure 8 contains similar plots showing the per- 
formance of the two Reed-Solomon codes used with 16-QAM against that of 
uncoded 16-QAM and uncoded QPSK. Figures 9 and 10 show the performance 
of the two Reed-Solomon codes used with both 32-QAM constellations against 
their respective uncoded cases. A summary of the BER performances of the error 
control codes over the corresponding uncoded cases is contained in Table 3 fol- 


lowing the performance graphs. 
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Figure 7. Uncoded Bit Error Rates: Top: BER vs. SNR. Bottom: BER vs. £,/\, 
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Table 3. ERROR CONTROL PERFORMANCE SUMMARY 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The simulation MFMLINK is the first implementation of error control cod- 
ing in multi-frequency modulation at NPS. As such, the codes utilized may not 
be the optimal codes for the given link pa meters, but do provide good insight 


into the benefits of error control coding for MR 





\ 


First, the enhanced error detection capabilities of the modified Reed-Solomon 


codes in the simulation support the introduction of some type of Automatic Re- 
peat Request (ARQ) strategy into the MFM environment. More specifically, a 
type 1 hybrid ARQ scheme is called for as described in Reference 2. Using this 
scheme with MFM, any correctable errors which are detected are corrected and 
if uncorrectable errors are detected, the receiver requests retransmission of par- 
ticular erroneous codewords or bauds. The proper combination of ARQ and FEC 
can provide a higher information throughput than one with ARQ alone, and 
provide higher system reliability than one with FEC alone [Ref. 2: p.478}. 
Second, although the decoding time for any given codeword or baud was not 
measured, the use of the decoding algorithms presented for the modified Reed- 
Solomon codes (which are solely based on the calculated syndromes) are simpler 
to implement in software than the iterative algorithms for conventional codes and 
felt to be less time consuming to run. The loss in information rate due to surren- 
dering an information symbol to parity is well worth the time saved in the de- 


coding process. 
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Finally, based on the graphs presented at the end of Chapter 4, the coding 
gains obtained at bit error rates of 10-> and 10-* for each code has proven the 
value of the codes used. Large gains were evident using error control with 
16-QAM modulation at a moderate cost in code rate. Smaller gains were evident 
using error control with 32-QAM modulation, particularly with the Gray-coded 
constellation, but at only a slight reduction in code rate. In any case, the loss in 
code rate due to the addition of parity symbols for error control is worth the 
coding gains achieved. 

Since this simulation was the first venture into MFM error control imple- 
mentation at NPS, several alterations to the simulation may provide improved 
error control performance. One might be to change the methods of channel de- 
coding, particularly in the case of 32-QAM. The approach used with 16-QAM 
may have improved performance over the maximum-likelihood approach. 
Constellation geometries with elements farther apart may improve performance, 
as seen in the results obtained with the 32-QAM constellations. Another alter- 
ation might be to vary the modulation technique used with the error control codes 
in the simulation. For example, using QPSK and the length 15 Reed-Solomon 
codes with each tone carrving two bits, two consecutive tones would represent one 
element of GF(2*) and 30 tones would represent a complete codeword. Interleav- 
ing the codeword symbols between bauds is another possible alteration, e.g., 
having a 16 tone information band with each tone representing one symbol from 
16 separate codewords would correspond to the transmission of 16 codewords in 
a 15 baud packet. Other methods of error control are also of interest, particularly 


the use of convolutional ’trellis codes and concatenated codes. 
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APPENDIX A. BINARY REPRESENTATIONS OF GALOIS FIELDS 


Table 4. ELEMENTS OF GF( 2*) 


Primitive Polynomial: 
P(x) =1+x4+ x4 
t 







Elemen Binary | 
Representation 






Table 5. ELEMENTS OF GF( 2) 


Primitive Polynomial: — 
P(x) = 14+ x? + x7 


Element Binary 
Representation 


ae ee 00000 
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APPENDIX B. SUPPORTING PROGRAMS 


A. GALOIS 


Crevededevedededededeve cede dete dete te reve Ne dee Fe Fe Fe Fe VEN TCI VETER TEI TE TE VET TE TVET ETE IE TE TE TENE TEPC ETE TENE FETE FEVER TENE TE FETE HOTCNE IE 
C THIS PROGRAM WILL READ THE M-BIT REPRESENTATIONS OF GF(2%*M) FROM * 


C A DATA FILE, THEN CALCULATE THE ADDITION TABLE ((N+1)*(N+1)) AND * 
C WRITE THE RESULTS INTO A DATA FILE GFA(Nt1). * 
Cc * 
Cievedede deve dene doce ie de dese se dete Tete ge Te Fe TENE TE TEN TENE TE EE ETE TE TEE NE ETE TENE TENE FETE TENE TE TE VE DEDEDE DEBE TENE ETE TEE TE TER IE TIE NE 


INTEGER P(0:15),G(0:15,0:7),A(0: 15,0: 15,0: 7),GFA(0: 15,0: 15) 
INTEGER S,N,M 


CALL EXCMS('FILEDEF 8 DISK GF16 DATA (PERM') 
OPEN( 7, FILE='GFA16' ,ACCESS=' SEQUENTIAL’ , FORM='FORMATTED' , STATUS= 


+ "UNKNOWN ) 
G OPEN(8,FILE='GF16' ,ACCESS='SEQUENTIAL' , FORM='FORMATTED' , STATUS 
C + ='UNKNOWN' ) 
Weis 
M=4 


C READ IN THE M-BIT REPRESENTATIONS FROM DATA FILE 


DO 10 I = 0,N 
READ(8,99) P(1),(G(1,J),J=0,M=1) 
WRITE(5,99) P(I),(G(1,J) ,J=0,M-1) 
10 CONTINUE 
99  FORMAT(1X,12,412) 


C ADD THE M-BIT ELEMENTS... 


DO 20 I=0,N 
DO 30 J=0,N 
DO 40 K=0,M-1 

S=G(1,K)+G(J,K) 
A(I,J,K)=MOD(S, 2) 

40 CONTINUE 

30 CONTINUE 

20 CONTINUE 


C NOW FIGURE OUT WHICH POWER OF ALPHA THESE M BITS REPRESENT... 


DO 50 J=0,N 
DO 60 J=0,N 
DO 70 K=0,N 
DO 80 L=0,M-1 
IF (ACI,J,L).NE. (G(K,L))) THEN 
GO TO 70 
ENDIF 





80 CONTINUE 
GFACI,J)=P(K) 
GO TO 60 

70 CONTINUE 

60 CONTINUE 

50 CONTINUE 


Cc NOW WRITE THESE POWERS OF ALPHA INTO THE ADDITION TABLE... 


DO 90 I=0,N 


WRITE(7,900) (GFA(I,J),J=0,15) 
THE FOLLOWING LINES ARE FOR LARGER FIELD CALCULATIONS... 

WRITE(7,900) (GFACI,J),J=16,31) 

WRITE(7,900) (GFA(1,J),J=32,47) 

WRITE(7,900) (GFA(1,J),J=48,63) 

WRITE(7,900) (GFA(I,J) ,J=64,79) 

WRITE(7,900) (GFA(I,J),J=80,95) 

WRITE(7,900) (GFA(I,J),J=96,111) 

WRITE(7,900) (GFA(I,J),J=112,127) 

WRITE(7,900) (GFA(I,J),J=128, 143) 

WRITE(7,900) (GFACI,J) ,J=144, 159) 

WRITE(7,900) (GFA(I,J) ,J=160,175) 

WRITE(7,900) (GFA(1,J),J=176,191) 

WRITE(7,900) (GFA(I,J) ,J=192,207) 

WRITE(7,900) (GFA(I,J),J=208,223) 

WRITE(7,900) (GFA(I,J),J=224,239) 

WRITE(7,900) (GFA(1,J),J=240,255) 
90 CONTINUE : 
900 FORMAT(16(1X,12)) 


AAAAANANRMAIaANCRaANAANANA 


ENDFILE( 7) 
CLOSE(7) 
END 








B. BITDIFF 
Cut ‘eve fot 
¢ THIS 'PROGRAN READS IN THE BIT REPRESENTATION OF A GALOIS FIELD AND * 


C 
C 
Cuvedst 


seve: eves! Seteveds ‘evs lowe se ves's alas! Pot es 'v23 Seve evedeses' eds vesteds' fe Seve eds eseves' ‘evevs estes! feve vet! teve see! ‘eveve ses fevey eves sede 


CALCULATES THE # OF BITS IN WHICH EACH ELEMENT DIFFERS. THEN IT * 
Riecoe ae pou a eo ic gos USE IN eidacen el ee Loe Cae ae 


seyevedsus seeeses'¢ stes' sles'e 


INTEGER ALPHA(0: 31) ,BIT(0: 31,0: 5),GFD(0: 31,0: 31),N,ROWS,DIFF 
CALL EXCMS('FILEDEF 9 DISK GF32 DATA (PERM' ) 
OPEN (7,FILE='GFD32' ,FORM='FORMATTED' , ACCESS=' SEQUENTIAL’ , 


+ STATUS=' UNKNOWN’ ) 


10 


40 


30 
20 


50 
1000 


N=31 
POWER=5 
ROWS=(N4+1)/16 


DO 10 I=0,N 
READ(9,*)ALPHA(1),(BIT(I,J) ,J=0, POWER-1) 
CONTINUE 
DO 20 I=0,N 
DO 30 J=0,N 
DIFF=0 
DO 40 K=0,POWER-1 
IF (BIT(I,K).NE.BIT(J,K)) THEN 
DIFF=DIFF+1 
ENDIF 
CONTINUE 
GFD(1,J)=DIFF 
CONTINUE 
CONTINUE 


DO 50 I=0,N 
WRITE(7,1000)(GFD(1,J) ,J=0,15) 
WRITE(7,1000)(GFD(1,J) ,J=16, 31) 

CONTINUE 

FORMAT(16(12,1X)) 

END 
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2 


TRACE 


Cwttet Ete TERE TT Te eT Te TS TET EE EEE ETT ETE ELE LEE TET EERE TOTES EP EPEDE PEELED LEE BEPEDEDED ELINA AAAS 
C THIS a ae ee FINDER PROGRAM * 
Cre? ‘ove sve Fedetedededetedededededetete face teee Seeded TOT Tete See Se CONC SESE EO 
C 

c { 


INTEGER A,B,C,D,GFA(0: 31,0: 31),N 


CALL EXCMS('FILEDEF 9 DISK GFA32 FILE (PERM' ) 
OPEN(6,FILE='TR32' , FORM='FORMATTED' , ACCESS='SEQUENTIAL' , 
+ STATUS=' UNKNOWN’ ) 


READ IN GF(2***M) ADDITION TABLE... 
N=31 


aga 


DO 10 I=0,31 
READ(9,*) (GFACI,J),J=0,15) 
READ(9,*) (GFACI,J) ,J=16,31) 
10 CONTINUE 


WRITE(6, 1001) 
1001 FORMAT('0',1X,'0') 
DO 55 I=2,31 
A=MOD( 2**(I-1),N)+1 es 
A=GFA(I,A) 
B=MOD(4*(I-1),N)+1 
B=GFA(A,B) 
C=MOD( 8**(I-1) ,N)+1 
C=GFA(B,C) 
D=MOD( 16*(I-1),N)+1 
WRITE(6,1002)GFA(C,D) 
55 CONTINUE 
1002 FORMAT(1X,12) 
END 


D. QAMCOM 


Cisietes Sevededete ede le tele te de Rede te Eee Pee S EEE DENELEEL IMENT ADIDAS 


C acne COMPUTER- THIS PROGRAM READS IN MAGNITUDE AND PHASE INFO ¥ 
C ON A QAM CONSTELLATION AND CONVERTS IT TO A REAL AND * 
C IMAGINARY PART WHICH IS STORED IN A DATA FILE "QANRI" ¥ 
C WHICH Is READ BY THE MEM LINK PROGRAM. ve 
Crretetiodedededede? Dedede terete te te te Sete Te Te Te TE SELES PETE EEE EEE LEE DEBSPE PEELE EEE DEDEDE EI ESEDEDE SEES. 


INTEGER ALPHA(0: 31) ,MAG(0: 31) ,N 
REAL RECO: 31), IM(0: 31) ,PI, PHASE(0: 31) 


CALL EXCMS('FILEDEF 9 DISK QAM32A DATA (PERM' ) 
OPEN(6,FILE='QAMRI32A' , FORM="FORMATTED' , ACCESS='SEQUENTIAL' , 
+ STATUS='UNKNOWN' ) 


N=31 
PI=4. O*ATAN( 1. 0) 


DO 1 I=0,N 
READ(9 ,**)ALPHA(1) ,MAG(I) , PHASE(I) 
PHASE(I)=PHASE(1)**PI/180. 0 
RE( 1 )=MAG(I)**COS( PHASE(I)) 

IM( I )=MAG(1)**SIN( PHASE(I)) 
WRITE(6,*)RE(1),1IM(I) ,MAG(I)/1. 0 
1 CONTINUE 

ENDFILE(6) 

CLOSE(6) 

END 
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APPENDIX C. SIMULATION PROGRAM: MFMLINK 


Credevedetede deve dese resets ied e Fete TeP Ne Nee Ne TENE TC FE ETE TE ETE TE ETE LELE TE TENE DE LEE DELETE LE EE TONED EE LEE DOLE TE GE PETE TEE PEER IC. 
eee ede dedevededede ve vedic de de Fete de Fete de Fede TEN METER TOT TOD TEETE TEE TE TE TCT TE TC T NE GE sededodedededsteketeddedieds 
Ce deve 
c* MFM COMMUNICATION LINK tere 
C* ads 
c* deve 
c* -THIS SIMULATES AN MFM COMMUNICATION LINK WHICH IS week 
c* CORRUPTED BY AWGN. IN DEALING WITH GALOIS FIELDS wee 
Cc* IT WAS NECESSARY TO TREAT THE ELEMENTS IN TERMS weve 
C* OF THEIR POWER OF ALPHA PLUS ONE, EXCEPT FOR THE weve 
c* ELEMENT "0" WHICH IS ALWAYS WRITTEN AS "0". ere 
C* seve 
Cirtevededevedededede deve de ve fede Fe eve Ne Te Ve NE ETL E TE TENE TE TE OLE LE TE TEE BEVELED TE TE TEE SEDGE PETE EE TELE DEDEDE TELE ELE DET E IETS 
c* weve 
C* -THE IMPORTANT ARRAYS/VARIABLES ARE AS FOLLOWS: ‘evs 
C* seve 
C* ARRAYS: D........... RANDOM DATA VECTOR (CODED CASES) vere 
Cc* DING oie ates RANDOM DATA VECTOR (UNCODED CASES) wee 
c* CWiisveie Sievecs @0s REED-SOLOMON CODEWORD FOR D Wve 
c* EM e sbeaneveueges ANALOG TRANSMISSION (OUTPUT OF INV-FFT) dere 
ce AWGN... 26... NOISE VECTOR weve 
c* BN, 2h ssech tere aise ANALOG RECEPTION (TK + AWGN, OUT OF FFT) ** 
c* Riss Stew pte. eres RECEIVED VECTOR wee 
C* DOU TP nas en RECEIVED VECTOR (UNCUDED CASES) vere 
C* Bisco ieot doh er gtens ESTIMATED ERROR IN R weve 
c* Vie Rvs aie sadsin ESTIMATED CODEWORD(R + E) were 
C* DHAT........ ESTIMATED DATA FROM R wie 
C* GFACI,J)....GALOIS FIELD SUM OF SYMBOLS weve 
c* GFD(I,J)....NUMBER OF BITS THAT 1 AND J DIFFER IN 

c* Se See eka wt RECEIVED VECTOR SYNDROME wee 
C* devs 
Cisdedederedkdevedsdevecede dete dete ve dese te Fee Fe Fe eT We ETE TE Te Se TET E TS TET TE ENE TE TE TE TENE TE TE TELE EE TEE EE TEE TENE TE EEOC EE ELON GE TE 


INTEGER B,C,K,BYTERR, ERRORS, FIXED, BDS , RUB, NLESSK, MAX 

INTEGER N,M, ITER, ERASE, DATERR,CCWW(17,0: 31) ,DIN(0: 255) ,DOUT(0: 255) 
INTEGER S(-2:2),GFA(0: 31,0: 31),R(0: 31),V(0: 31) ,E(0: 31) 
INTEGER TRACE(0: 31) ,T4(0: 31) ,D(0: 31) ,CW(0: 31) ,DD(17,0: 31) 
INTEGER GAMMA1,GAMMA2,GAMMA3 , ZEROS , TRIERR, DUOERR , UNOERR , OOPS 
INTEGER DHAT(0: 31),BDTYPE,CODE,K1,K2,PTS,CPB,PTR, BAUDNO 
INTEGER GFD32(0: 31,0: 31) ,DBERR,CWBERR,GFD16(0: 15,0: 15) 
INTEGER QPDIFF(0: 3,0: 3) 

REAL D1(0: 31),DO(0: 255) ,QaMRE(0: 31) , QAMIM(0: 31) , QANMAG( 0: 31) 
REAL MEAN,NVAR, AWGN( 0: 4095) , TEMP1(0: 4095) , TEMP2(0: 4095) 
COMPLEX QAMCOM(0: 31),TX(0: 4095) ,RX(0: 4095) ,T(0: 4095) 

REAL SUM,SUMSQ,SNR,NFIX,QPRE(0: 3) ,QPIM(0: 3) 


CALL EXCMS('FILEDEF 9 DISK GFA16 DATA (PERN') 
CALL EXCMS('FILEDEF 7 DISK TR16 DATA (PERM)"') 
CALL EXCMS('FILEDEF 8 DISK QAMRI16 FILE (PERM' ) 
CALL EXCMS('FILEDEF 1 DISK GFD16 FILE (PERM') 
CALL EXCMS('FILEDEF 11 DISK GFA32 FILE (PERM') 





C SELECT WHICH CODE DESIRED... 


413 


21 





CALL EXCMS('FILEDEF 12 DISK TR32 FILE (PERM)') 
CALL EMCMS('FILEDEF 13 DISK QAMRI32A FILE (PERM') 
CALL EXCMS('FILEDEF 31 DISK GFD32 FILE ( PERM’) 
CALL EXCMS('FILEDEF 14 DISK QPSKRI FILE (PERM') 
CALL EXCMS('FILEDEF 34 DISK QPSKDI FILE (PERM') 


OPEN(6,FILE='RSRSLTS' ,FORM='FORMATTED' ,ACCESS='SEQUENTIAL', 


+ STATUS=' UNKNOWN" ) 


OPEN( 33,FILE='C2C1NEW' , FORM='FORMATTED' , ACCESS='SEQUENTIAL' 


+ STATUS=' UNKNOWN’ ) 


WRITE(5, 


*) ' 
WRITE(5,*)'1. SELECT COD 
*) : 


WRITE(5, 
WRITE(5,*)' 
WRITE(5,*)' 
WRITE(5,7")' 
WRITE(5,%*)’ 
WRITE(5,*)' 
WRITE(5,*)! 
READ( 5, %*) CODE 


Auf wnr oO 


TYPE acs 


UNCODED 
(75,10) 
(31,26) 
(15,12) 
(31,28) 
UNCODED 
UNCODED 


(16-QAM)' 

REED-SOLOMON CODE (16-QAM)' 
REED SOLOMON CODE (32-QAM)' 
REED SOLOMON CODE (16-QAM)' 
REED SOLOMON CODE (32-QAM)' 
(32-QAM)' 

QPSk' 


IF ((CcDE. LT. 0). OR. (CODE. GT. 6)) THEN 
WRITE(5,%*)'BAD CHOICE OF CODE--TRY AGAIN?' 


GO TO 413 
ENDIF 
IF (CODE. EQ.6) THEN 
WRITE(5,*)' ' 


WRITE(5,*)'2. SELECT BAUD TYPE (1-5)' 
E 


READ(5 .*) BDTYP 


IF ((BDTYPE. GT. 5).OR. (BDTYPE. LT. 1)) THEN 


GO TO 21 
ENDIF 
ENDIF 


IF (( CODE. EQ. 1). OR. (CODE. EQ. 0). OR. (CODE. EQ. 3)) THEN 


N=15 

MAN=15 

NLESSK=5 

IF (CODE.EQ.3) THEN 
MAN=15 
NLESSK=3 

ENDIF 

WRITE(5,*)' 


WRITE(5,**)'2. SELECT BAUD TYPE (1-5)' 


READ(5 ,**) BDTYPE 


IF ((BDTYPE. GT.5).OR. (BDTYPE. LT. 1)) THEN 


GO TO 2 
ENDIF 
ENDIF 
IF (CODE. EQ.6) THEN 
MAX=1 
ENDIF 











IF ((CODE. EQ. 2). OR. (CODE. EQ. 4). OR. (CODE. EQ.5)) THEN 
N=31 
MAX=31 
NLESSK=5 
IF (CODE. EQ.4) THEN 
MAX=31 
NLESSK=3 
ENDIF 
WRITE(5,*)' ' 
221 WRITE(5,%*)'2. SELECT BAUD TYPE (2-5)' 
READ(5 ,*)BDTYPE 
IF ((BDTYPE. GT. 5). OR. (BDTYPE. LT. 2)) THEN 


















GO TO 221 

ENDIF 

ENDIF 

IF ((CODE. GT. 0). AND. (CODE. LT.5)) THEN 
WRITE(5,*)' 
WRITE(5,*)'3. ERASE CODEWORDS WITE UNCORRECTED ERRORS es 
WRITE(5,*)' ENTER "0" TO LEAVE CODEWORDS AS IS' 
WRITE(5,*)' ENTER "1" TO ERASE’ 
READ(5 ,**)RUR 

ENDIF 






WRITE(S ,*)) / 
WRITE(5,*)'4. SELECT SNR(DB,REAL NUMBER)' 
READ(5,**)SNR 





aAANQ 







C READ IN DATA AND CONSTANTS... 
C 






CALL BAUDS(CODE,BDTYPE,K1,K2,PTS,CPB,M) 
IF (CODE. EQ.6) THEN 
DO 108 1=0,3 
READ( 34, **)(QPDIFF(I,J) ,J=0, 3) 
WRITE(5 ,**)(QPDIFF(T,J) ,J=0,3) 
READ( 14,*)QPRE(1) ,QPIM(1) 
QAMCOM( I) =CMPLX( QPRE( I) ,QPIM(I)) 












108 CONTINUE 
GO TO 597 
ENDIF 
CALL DATA(CODE,N,GFA, TRACE ,T4,QAMCOM , QAMRE , QAMIM, QAMMAG ,GFD16, 
+ GFD32) 








597 NFIX=10*LOG10( (PTS/2)/(K2-K1)/1.0) 







DO 544 SNR=15.0,22.0,1.0 





BYTERR=0 
NOERR=0 
UNOCERR=0 
ICOERR=0 
‘I RIERR=0 
FIXED=0 
DATERR=0 















CWBERR=0 
DBERR=0 
ITER=0 
ERASE=0 
MEAN=0. 0 


DO 543 BDS=1,100 


DO 15 I=0,MAX-1 
CW(1)=0 
D(1I)=0 
15 CONTINUE 
DO 16 I=0,PTS-1 
TX(1)=(0. 0,0. 0) 
RX(I)=(0. 0,0. 0) 
T(I)=(0. 0,0. 0) 
16 CONTINUE 


C MAKE UP RANDOM DATA, ENCODE IT, LOAD IT IN TRANSMITTER... 
IF ((CODE. EQ. 0). OR. (CODE. EQ.5)) THEN 
CALL NOCODE(CPB,DO,T,K1,PTS,QAMCOM,DIN,N) 
GO TO 797 
ENDIF 
IF (CODE. EQ. 6) THEN 
CALL QPSK(DO,K1,K2,PTS,QAMCOM,DIN,T) 
GO TO 797 
ENDIF 
PTR=K1 
DO 776 J=1.CPB 
CALL RNUN( MAX-NLESSK,D1) 
DO 2 I=0,MAX-NLESSK-1 
D1(1)=INT( (N#1)*D1(1)) 
1 CONTINUE 
CALL INTRAN(D1,D,MAX,NLESSK) 
CALL ENCODE(N,CW,D,GFA, CODE ,NLESSK , MAX) 
DO 17 I=PTR,PTR+MAX-1 
T(1)=QAMCOM( CW(I-PTR)) 
T( PTS-I)=CONJG(T(1)) 
DD( J, 1-PTR)=D( 1-PTR) 
CCWW(J,1-PTR)=CW(I-PTR) 
17 CONTINUE 
PTR=PTR+MAX 
776 CONTINUE 


C TRANSMIT THE CODEWORD... 
797 CALL BITREV(PTS ,M,T,TX) 
CALL INVFFT( PTS,M,TX) 


C ANY NOISE IN COMMUNICATION CHANNEL? (R(I) IS RECEIVED VECTOR)... 
SUMSQ=0. ) 
SUM=0. 0° 
DO 13 1=0,PTS-1 





TXCI)=TXC1) /2+CONJG(TX(1I))/2 
SUM=SUM+REAL(TX( 1) ) 
SUMSQ=SUMSQt+( REAL( TX( I) )**2) 
13 CONTINUE 
NVAR=SUMSQ/( PTS-1) /( 10**( (SNR-NFIX)/10.0)) 


CALL GAUSS(MEAN,NVAR, PTS ,AWGN,TEMP1, TEMP2) 
DO 617 [=0,PTS-1 
TXCI)=TXCI)+AWGNCT) 
617 CONTINUE 


C RECEIVE THE DATA... 
CALL BITREV( PTS ,M,TX,RX) 
CALL FFT(PTS,M,RxX) 
C DECODE EACH RECEIVED CODEWORD SEQUENTIALLY... 
IF (CODE. EQ.0) THEN 
CALL WHAT16(RX,CPB,PTS,K1,DOUT) 
PTR=K1 
GO TO 667 
ENDIF 
IF (CODE.EQ.5) THEN 
CALL what32(RX,N, PTS, QAMRE,QAMIM,DOUT,K1,CPB) 
PTR=K1 
GO TO 668 
ENDIF 
IF (CODE. EQ. 6) ‘THEN 
CALL WHATQP(RX,K1],K2,PTS,DOUT) 
GO TO 669 
ENDIF 


PTR=K1 
DO 777 J=1,CPB 
B=0 
C=0 
K=0 
DO 881 I=0,MAX-1 
CW(1)=CCKW(J,1) 
D(1)=DD(J,1) 
DHAT(I)=0 
V(I)=0 
E(I)=0 
R(1)=0 
881 CONTINUE 
IF ((CODE. EQ. 1). OR. (CODE. EQ. 3)) THEN 
CALL WHAT16(RX,MAX,PTS,PTR,R) 
ENDIF 
IF ((CODE. EQ. 2). OR. (CODE. EQ.4)) THEN 
CALL what32(RX,MAX, PTS , QAMRE, QAMIM,R, PTR, MAX) 
ENDIF 
PTR=PTR+MAX 


C CALCULATE THE SYNDROME... 
ZEROS=0 
NEAKEN=0 





OOPS=0 
CALL SYNDRO(R,S,N,NEAREN,OOPS ,GFA, ZEROS , MAX, CODE) 


Cc CHECK FOR NO ERRORS... 
199 IF (ZEROS. EQ. NLESSK) THEN 
NOERR=NOERR+1 { 
CALL RUBOUT(RUB,R,V,DHAT,N,NLESSK, ZEROS , MAX) 
GO TO 112 


c GET INTO REST OF DECODING ALGORITHM... 
IF ( (CODE. EQ. 1). OR. (CODE. EQ. 2)) THEN 
BAUDNO=J 
B=0 
c=0 
kK=0 
CALL DECTWO(R,S,V,DHAT,E,ERRORS , UNOERR , DUOERR,TRIERR 
+ »FIXED,B,C,K,ERASE,TRACE,N,GFA,GAMMA1 , GAMMA2 , GAMMNA3 
+ » RUB, NLESSK, ZEROS , CODE , MAX) 
ENDIF 
IF (( CODE. EQ. 3). OR. (CODE. EQ. 4)) THEN 
CALL DECONE(S,V, ZEROS ,RUB,R,DHAT,N,NLESSK,MAX,GFA,CODE, 
+ UNOERR , DUOERR) 


C NOW CHECK IF ANY SYMBOL ERRORS WERE MADE (COMPARE ESTIMATE 
Cc OF CODEWORD SENT TO ACTUAL CODEWORD SENT)... 
112 ERRORS=0 
BO 61 1=0,MAX-1 
IF (VC1I).NE.CWCI)) THEN 
BYTERR=BYTERR+1 
ERRORS=1 
IF ((V(I). EQ. -1). AND. ((CODE. EQ. 1). OR. (CODE. EQ. 3))) 
+ THEN 
CWBERR=MAXN*4 
DBERR=(MAX-NLESSK)*4 
SO TO 388 
ENDIF 
IF ((V(I). EQ. -1). AND. ((CODE. EQ. 2). OR. (CODE. EQ. 4))) 
+ THIN 
CWBERR=MAX*5 
DBERR=(MAX-NLESSK)*5 
GO TO 388 
ENDIF 
IF ((CODE. EQ. 0). OR. (CODE. EQ. 1). OR. (CODE. EQ. 3)) THEN 
CWBERR=CWBERR+GFD16(V(1I).CW(1)) 
GO TO 133 
ENDIF 
CWBERR=CWBERR+GFD32(V(1) ,CW(1)) 
133 ENDIr 
ERROTS=1 
61 CONTINUE 
DO 336 1=0,MAX-NLESSK-1 


49 








IF (D(1I).NE. DHAT(I)) THEN 
DATERR=DATERR+1 
IF (DHAT(I).EQ.-1) THEN 
DBERR=DBERR+4 
ELSE 
IF (( CODE. EC. 0). OR. (CODE. EQ. 1). OR. (CODE. EQ. 3)) THEN 
DBERR=DBERR+GFD16(D(1) ,DHAT(1)) 


. GO TO 231 
ENDIF 
DBERR=DBERR+GFD32(D(1) ,DHAT(I)) 
231 ENDIF - 
ENDIF 
336 CONTINUE 
Cc fa ae ak as Stab sma “ess aa as as 0 0st a laa tgs ts ss a cnt mapa at ew aa ew eve tae ie sem js see es es ee ae ie 
388 BAUDNO=J 
C IF (ERRORS. EQ. 1) THEN 
C113 CALL OUTPUT(S,GAMMA1,GAMMA2 ,GAMMA3,R,V,B,C,K,CW,D, 
Ge BAUDNO , DHAT, BDTYPE , MAX) 
C ENDIF 
335 NOP=0 
99 TER=ITER+1 
777 CONTINUE 
667 IF (CODE.EQ.0) THEN 


DO 942 I=0,CPB-1 
IF (DIN(I).NE. DOUT(I)) THEN 
DATERR=DATERR+1 
DBERR=DBERR+GFD16(DIN( I) ,DOUT(I)) 


ENDIF 
942 CONTINUE 
ENDIF 
668 IF (CODE. EQ.5) THEN 


DO 944 1=0,CPB-1 
IF (DIN(1).NE.DOUT(1)) THEN 
DATERR=DATERR+1 
DBERR=DBERR+GFD32(DIN(1) ,DOUT(I)) 


ENDIF 
944 CONTINUE 
ENDIF 
669 IF (CODE. EQ.6) THEN 


BO 945 I=0,K2-K1 
IF (DINCI).NE. DOUT(I)) THEN 
DATERR=DATERR+1 
DBERR=DBERR+QPDIFF(DIN(I) ,DOUT(TI)) 
ENDIF 
945 CONTINUE ‘ 
ENDIF . 
543 CONTINUE 
IF (CODE. EQ. 0) THEN 
BE=CPB*4*(BDS-1)/1.0 
BE=DBERR/BE 
WRITE(33,*)SNR, BE 
ENDIF 
IF (CODE. EQ.5) THEN 
BE=CPB**5*(BDS-1)/1.0 
BE=DBERR/BE 
WRITE(33,*)SNR,BE 








544 


ANMANAAANRAANAAARAAAAAAANAAAAAAAAARANAAANANANANNAAN 
“sl 
ur 
a 





ENDIF 

IF (CODE. EQ. 6) THEN 
BE=(K2-K1+1)*2%(BDS-1)/1.0 
BE=DBERR/BE 
WRITE( 33 ,**) SNR, BE 

ENDIF 

IF (( CODE. EQ. 1). OR. (CODE. EQ. 3)) THEN 
BE=CPB*N*4*( BDS-1)/1.0 
BE=DBERR/BE 
WRITE(33,*)SNR,BE 

ENDIF 

IF ( (CODE. EQ. 2). OR. (CODE. EQ.4)) THEN 
BE=CPB*N*5*(BDS-1)/1.0 
BE=DBERR/BE 
WRITE(33,*)SNR,BE 


ENDIF 

CONTINUE 

WRITE(6,*)' SNR (DB) =',SNR 
WRITE(6,7*)' BAUD TYPE = — ,BDTYPE 


IF (CODE. GT.0) THEN 

IF (CODE. EQ. 1) THEN 
WRITE(6,*)' CODING CHOICE 
GO TO 756 

ENDIF 

IF (CODE. EQ. 2) THEN 
WRITE(6,*)'CODING CHOICE = REED-SOLOMON (31,26)' 
GO TO 756 

ENDIF 

IF (CODE. EQ. 3) THEN 
WRITE(6,**) ‘CODING CHOICE = REED-SOLOMON (15,12)' 
GO TO 756 

ENDIF 

IF (CODE. EQ.4) THEN 
WKLTE(6,%*)'CODING CHOICE = REED-SOLUMON (31,28)' 


REED-SOLOMON (15,10)' 


GO TO 756 
ENDIF 
WRITE(6,*)' CW SYMBOL ERRORS AFTER CORRECTION=' ,BYTERR 
WRITE(6,*)' TOTAL SYMBOLS TRANSMITTED=' ,N* ITER 
WRITE(6,*)' RESULTING CODEWORD BIT ERRORS=' ,CWBERR 
IF ((CODE. EQ. 1). OR. (CODE. EQ. 3)) THEN 

WRITE(6,*)' CODEWORD BITS TRANSMITTED=' ,N*ITER*4 


BE=N*\ITER*4/1. 0 
BE=CWBERR/ BE 
ELSE 
WRITE(6,**)' CODEWORD BITS TRANSMITTED=" ,N*ITER*5 
BE=N*ITER*5/1. 0 
BE=CWBERR/BE 


ENDIF 
WRITE(6,*)' CODEWORDS WITH NO ERRORS=' ,NOERR 
WRITE(6,*)' DETECTED SINGLE ERRORS=' ,UNOERR 
WRITE(6,*)' DETECTED DOUBLE ERRORS="' ,DUOERR 


IF (( CODE. EQ. 1). OR. (CODE. EQ. 2)) THEN 
WRITE(6,*)'’ DETECTED TRIPLE (OR MORE) ERRORS=',TRIERR 
WRITE(6,**)' CORRECTED TRIPLE (BURST) ERRORS=' , FIXED 


St 


AQAAAANANAAAAAANAANANANAANANAN 


AQAaAGAIAANAAANA 


he 


slealeclectee's 
vevevedks's 


ENDIF 
WRITE(€,%*)' ERASURES=' , ERASE 
WRITE(6 ,**)' TOTAL CODEWORDS TRANSMITTED=' ,ITER 
ENDIF 
IF (CODE. EQ.0) THEN 
WRITE(6,*) ‘CODING CHOICE = UNCODED (QAM-16)' 


WRITE(6,%*)' DATA SYMBOL ERRORS=' ,DATERR 

WRITE(6,**)' DATA SYMBOLS TRANSMITTED=' ,CPB*( BDS-1) 

WRITE(6,**)' DATA BIT ERRORS=' ,DBERR 

WRITE(6,*)' DATA BITS TRANSMITTED="' ,CPB*4**(BDS-1) 
ENDIF 


IF (CODE.EQ.5) THEN 
WRITE(6,%) ‘CODING CHOICE = UNCODED (QAM-32)' 


WRITE(6,*)' DATA SYMBOL ERRORS=' ,DATERR 
WRITE(6,*)' DATA SYMBOLS TRANSMITTED=' ,CPB*( BDS-1) 
WRITE(6,*)' DATA BIT ERRORS=' ,DBERR 
WRITE(6,**) ' DATA BITS TRANSMITTED=' ,CPB*5**(BDS-1) 

ENDIF 

WRITE(6,*)' ' 

WRITE(6,*)'BIT ERROR RATE = ',BE 

WRITE(6,*)' TOTAL BAUDS EXAMINED=' ",BDS-1 


WRITE(33,**)' TOTAL BITS IN EACH SNR=' ,(K2-K1)*4**(BDS-1) 
ENDFILE( 6) 

CLOSE(6) 

ENDFILE(33) 

CLOSE( 33) 

END 


tevededevevetede sete secede Tose cece Tete ce TOTNES NCTE TESTE VTE EN RESERVED ENCE SES ER SEE SNE TET EDCLES ES CNC TES 


SUBROUTINE BAUDS(CODE,BDTYPE,K1,K2,PTS,CPB,M) 


SUBROUTINE TO SET SOME INITIAL VALUES, DEPENDING ON BAUD TYPE... 


CONSTANTS: PTS..... FFT SIZE 
js Oaoerererea POWER OF 2 (OF PTS) 
Kliaucce’ FIRST TONE NUMBER FOR 16-20KHZ BAND WITH 
SAMPLING FREQ=61440HZ AND FFT SIZE PTS 
| erate LAST TONE FOR THIS BAND 
CPB..... CODEWORDS PER BAUD (ONE (15,10) REED-SOLOMON 


CODEWORD PER 15 TONES IN THIS BAND) 
INTEGER BDTYPE,K1,K2,PTS,CPB,M,CODE 


IF (BDTYPE. EQ. 1) THEN 
PTS=256 
M=8 
CPB=1 
K1=68 
K2=83 
ENDIF 
IF (BDTYPE. EQ. 2) THEN 
PTS=512 
M=9 
CPB=2 


mn 
w 





K1=135 
K2=166 
ENDIF 
IF (BDTYPE.EQ.3) THEN 
PTS=1024 
M=10 
CPB=4 
K1=269 
K2=332 
ENDIF 
IF (BDTYPE.EQ.4) THEN 
PTS=2048 
M=11 
CPB=8 
K1=537 
.2=664 
ENDIF 
IF (BDTYPE.EQ.5) THEN 
PTS=4096 
M=12 
IF (CODE. EQ. 1) THEN 
CPB=17 
ELSE 
CPB=16 
ENDIF 
K1=1073 
K2=1328 
ENDIF 
IF ( (CODE. EQ. 0). Oh. (CODE. GE.5)) THEN 
CPB=K2-Ki+1 
ENDIF 
IF (( CODE. EQ. 2). OR. (CODE. EQ.4)) THEN 
CPB=INT(CPB/2) 
ENDIF 
RETURN 
END 


Crededededede dette dete de de terete te dete eS TE BETTE NE TEE TETE TENET TEENIE NE TE TET ETE RENE TE TEE DELETE DEDEDE NE TEE DENTE ETE TOTES 
SURROUTTNE RS Paes “GFA TRACE. Ti. TON GM Gin On tie. 
+ GFD16,GFD32) 


C THIS SUBROUTINE PROVIDES DATA TO THE PROGRAM DEPENDING ON THE CHOICE 
C OF CODE. 

Cc 

Cc ROWS. ..... NUMBER OF ROWS IN DATA FILE FOR EACH ROW IN GF 

C ADDITION TABLE (16 ENTRIES IN EACH ROW OF DATA FILE) 
C QAMRE..... REAL PART OF A QAM CONSTELLATION ELEMENT 

C QAMIM..... IMAGINARY PART OF A QAM CONSTELLATION ELEMENT 

C QAMCOM.... THE COMPLEX ARRAY OF (QAMRE, QA'MRI) 

C QAMMAG....MAGNITUDE OF THAT QAM ELEMENT 

Cc 

Cs aia ein Se Deen cain ae Vere ws ee ewe tees eee ce ee bee etctoccl ect ceceees 


INTEGER N,SUM, ROWS, CODE 

INTEGER GFAC 0: N,0:; N), GFD16(0: 15,0: 15) ,TRACE( 0: N) ,T4(0: N) 
INTEGER GFD32(0: 31,0: 31) 

REAL QAMRE( 0: N) ,QAMIM(0:N) ,X,¥Y,QAMMAG(0: N) 





COMPLEX QAMCOM(0:N) 


IF (CODE. EQ. 1) THEN 
READ(7 ,*)(TRACE(1I) ,I=0,N) 
READ(7 ,*)(T4(1) , I=0,N) 
ENDIF 
IF (CODE. EQ. 2) THEN 
READ( 12,*)(TRACE(1) ,I=0,N) 
ENDIF 
ROWS=INT(N+1)/16 
DO 10 I=0,N 
SUM=0 
DO 20 J=1,ROWS 
IF ((CODE. EQ. 0). OR. (CODE. EQ. 1). OR. (CODE. EQ. 3)) THEN 
READ(9 ,*)(GFAC I,K) ,K=0,15) 
READ(1,*)(GFD16(1,K),K=0,15) 
GO TO 20 
ENDIF 
READ(11,*) (GFACI,K) ,K=SUM,SUM+15) 
READ(31,*) (GFD32(1,K),K=SUM,SUM+15) 
SUM=SUM+16 


20 CONTINUE 
10 CONTINUE 
DO 40 I=0,N 


IF (( CODE. EQ. 0). OR. (CODE. EQ. 1). OR. (CODE. EQ. 3)) THEN 
READ(8 ,*)QAMRE( 1) ,QAMIM(I) , QAMMAG( I) 
GO TO 35 
ENDIF 
READ( 13 ,*)QAMRE( I) ,QAMIM(I) ,QANMAG( I) 
35 X=QAMRE( 1) 
Y=QAMIM(T) 
QAMCOM( I)=CMPLX(X, Y) 
40 CONTINUE 
RETURN 
BNP 


Cisdederedevedede cede vevede deve ve reese de sede rete ve de se Le TEFEN TE TENE EVE TENE NE TENCE TEE LENE LENE ELE VE TELE DEVE NE NE TELE TE TE VE PETE TOC TE TE 


SUBROUTINE NOCODE( CPB, D0, T,K1,PTS,QAMCOM,DIN,N) 


THIS SUBROUTINE GIVES RANDOM DATA FOR THE UNCODED CASES AND 
LOADS IT INTO THE TRANSMITTER. . 


DO........ REAL NUMBER EQUAL TO INTEGER DIN (FOR INV-FFT) 


C 

Cc 

Cc 

C DIN....... INTEGER REPRESENTING AN ELEMENT OF GF(2***M) 
Cc 

C Diack aided say eriene CONJUGATE SYMMETRIC SPECTRUM FOR INV-FFT 

C 


INTEGER CPB,K1,DIN(0: 255) ,PTS,N 
REAL DOC 0: 255) 
COMPLEX T(0: 4095) ,QAMCOM(0: 31) 


CALL RNUN(CPB ,DO) 
DO 10 I=0,CPB-1 
DOC I j)=INT( (N+1)*D0(1)) 
DO 20 J=0,N 
IF (DO(I).EQ.J/1.0) THEN 





DIN(I)=J3 
ENDIF 
20 CONTINUE 
T(K1+I)=QAMCOM(DIN(I)) 
T(PTS-K1-I)=CONJG(T(K1+I)) 
10 CONTINUE 
RETURN 
END 


Crevevedeverevedede tate de rete dete Lede eB F BETTE TEE ETE TELE L LENE ELLE TERETE BEDE LDPE NEDE DELLE LER LLL BELLE TETAS 


SUBROUTINE QPSK(DO,K1,K2,PTS,QAMCOM,DIN,T) 


THIS SUBROUTINE GIVES RANDOM DATA FOR UNCODED QPSK AND 
LOADS IT INTO THE TRANSMITTER... 


DO........ REAL NUMBER EQUAL TO INTEGER DIN (FOR INV-FFT) 


Cc 

Cc 

Cc 

Cc DING 8 a2.6:4% INTEGER REPRESENTING AN ELEMENT OF GF(2**M) 
Cc 

Cc Ts te eiencin ahs CONJUGATE SYMMETRIC SPECTRUM FOR INV-FFT 

Cc 


INTEGER K1,K2,DIN(0: 255) ,PTS 
REAL DOC0: 255) 
COMPLEX T(0: 4095) ,QAMCOM(0: 31) 


CALL RNUN(K2-K1+1.D0) 
DO 10 I=0,K2-K1 
DOC I)=INT(4*D0(1)) 
DO 20 J=0,3 
IF (POCI).EQ.J/1.0) THEN 
DINCI)=J 
ENDIF 
20 CONTINUE 
T(Ki+I)=QAMCOM(DIN(1)) 
T(PTS-K1-1)=CONJG(T(K1+I)) 
10 CONTINUE 
RETURN 
END 


Critededetorededevedededededededevede dered ve Petevedede de deve dete Tee Te Me TeV Te Le Te ENE Te TENE TTC NE NEN TE TERE TENE IER Te Te TE ee 


SUBROUTINE ENCODECN, cM, D,GFA,CODE, NLESSK. SMAX) 


C THIS SUBROUTINE IS THE REED-SOLOMON ENCODER (SYSTEMATIC FORM). 
C 

C GO-G4........ COEFFICIENTS OF GENERATOR POLYNOMIAL 

C SR ie eaee se SHIFT REGISTERS 

C 

C 


INTEGER N, INPUT ,GO,G1,G2,G3,G4,CODE ,NLESSK , MAX 
INTEGER CW(0: MAX-1) ,D(0: MAX-1) ,GFACO:N,0:N),SR(5) 


C SET GENERATOR POLYNOMIAL CONSTANTS... 
G1=5 


nN 
‘aA 





G2=12 
G3=12 
G4=5 
IF (CODE.EQ. 2) THEN 
G1=14 
G2=18 
G3=18 { 
G4=14 
ENDIF 
IF (CODE. EQ. 3) THEN 
GO=4 
G1=12 
G2=11 
ENDIF 
IF (CODE. EQ.4) THEN 
GO=4 
G1=13 
G2=12 
ENDIF 
C ENCODE THE DATA FOR THE SYSTEMATIC (N,N-5) CODE... 


411 DO 8 J=1,NLESSK 
SR(J)=0 
8 CONTINUE 
DO 9 J=0,MAX-NLESSK-1 
CWC MAX-1-J)=D( MAX-NLESSK-1-J) a 
INPUT=GFA(SR(1),D(MAX-NLESSK-1-J)) 
IF (INPUT. EQ. 0) THEN 
DO 40 K=1,NLESSK-1 
SRC K)=SR(K+1) 
40 CONTINUE 
SRCNLESSK)=INPUT 
GO TO 9 
ENDIF 
IF (NLESSK. EQ. 5) THEN 
SR(1)=GFACSR(2) ,MODC INPUT+G4-2 ,N)+1) 
SR(2)=GFACSR(3) ,MOD( INPUT+G3-2,N)+1) 
SRC 3)=GFACSR(4) ,MODC INPUT+G2-2,N)+1) 
SR(4)=GFACSR(5) ,MOD(C INPUT+G1-2,N)+1) 
SR(5)=INPUT 
ENDIF 
IF (NLESSK. EQ. 3) THEN 
SR(1)=GFA(SR(2) ,MOD(C INPUT+G2-2 ,N)+1) 
SR(2)=GFA(SR(3) ,MOD( INPUT+G1-2,N)+1) 
SR(3)=MOD( INPUT+G0-2,N)+1 
ENDIF 
9 CONTINUE 
DO 11 I=0,NLESSK-1 
CWC I)=SRCNLESSK-1) 
11 CONTINUE 
RETURN 
ae 
Cardedevedevede dace vedere dete seve terete tote te Be eT Bede Novedededededededededevede de tele te sede te Lege FORESEES LTP REPRE EEE, 
SUBROUTINE WHATIO(S, MAX, ‘L, PTR 2 R) 


56 








SUBROUTINE TO DECIDE WHICH ELEMENT OF GF(N+1) EACH ELEMENT 
OF THE RECEIVED DATA CORRESPONDS TO. 
PTR...... CURRENT CODEWORD STARTING TONE 
a4) arerereree AN ELEMENT FROM THE RECEIVER (FFT) 
MAX Nici che CODEWORD BLOCK LENGTH 
Liv.bs 5 o8%e.8%s FFT SIZE 
15 @) errr ee THE RECEIVED CODEWORD IN TERMS OF ELEMENTS OF GF(N+1) 


INTEGER MAX,R(0: MAX-1),PTR 
COMPLEX X(0: L-1) 
REAL A,B,2Z,MID 


MID=150. 0 
DO 10 J=PTR,PTR+MAX-1 
A=REAL(X(J)) 
B=AIMAG(X(J)) 
Z=CABS(X(J)) 
IF ((A.GE.0). AND. (B.GE.0)) THEN 
IF ((Z.GE.MID). AND. (A.GE.B)) THEN 
R( J-PTR)=4 
GO TO 10 
ENDIF 
IF ((Z.GE.MID). AND. (A. LT. B)) THEN 
R( J-PTR)=7 
GO TO 10 
ENDIF 
IF ((Z. LT. MID). AND. (A. GE. B)) THEN 
R( J-PTR)=0 
GO TO 10 
ENDIF 
R( J-PTR)=3 
GO TO 10 
ENDIF 
IF ((A.GE.0). AND. (B. LT.0)) THEN 
IF ((Z. GE. MID). AND. (A. GE. ABS(B))) THEN 
R( J-PTR)=15 
GO TO 10 
ENDIF 
IF ((Z.GE.MID). AND. (A. LT. ABS(B))) THEN 
R( J-PTR)=14 
GO TO 10 
ENDIF 
IF ((Z. LT. MID). AND. (A. GE. ABS(B))) THEN 
R( J-PTR)=1 
GO TO 10 
ENDIF 
R( J-PTR)=9 
GO TO 10 
ENDIF 
IF ((A. LT. 0). AND. (B. LT. 0)) THEN 
IF ((Z.GE. MID). AND. (ABS(A). GE. ABS(B))) THEN 
R(J-PTR)=8 
GO TO 10 





ENDIF 
IF ((Z.GE.MID). AND. (ABS(A). LT. ABS(B))) THEN 
R( J-PTR)=13 
GO TO 10 
ENDIF 
IF ((Z. LT. MID). AND. (ABS(A). GE. ABS(B))) THEN 
R( J-PTR)=5 : 
GO TO 10 
ENDIF 
R( J-PTR)=11 : 
GO TO 10 
ENDIF 
IF ((A.LT. 0). AND. (B.GE.0)) THEN 
IF ((Z.GE.MID). AND. (ABS(A).GE.B)) THEN 
R( J-PTR)=10 
GO TO 10 
ENDIF 
IF ((Z.GE.MID). AND. (ABS(A). LT. B)) THEN 
R( J-PTR)=12 
GO TO 10 
ENDIF 
IF ((Z. LT. MID). AND. (ABS(A).GE.B)) THEN 
R( J-PTR)=2 
GO TO 10 
ENDIF 
R( J-PTR)=6 
ENDIF 
10 CONTINUE 
RETURN 
END 


Seateate niente: 


SUBROUTINE WHATQP(X,K1,K2,L,R) 


C SUBROUTINE TO DECIDE WHICH ELEMENT OF GF(N+1) EACH ELEMENT 

C OF THE RECEIVED DATA CORRESPONDS TO. 

C 

Cc RG) ice yore AN ELEMENT FROM THE RECEIVER (FFT) 

C 19 FIRST TONE IN BAUD 

C | ae SECOND TONE IN BAUD 

C De iiacie eae FFT SIZE 

C RO) ae bass THE RECEIVED CODEWORD IN TERMS OF ELEMENTS OF GF(N+1) 
C 


INTEGER RCO: K2-K1+1) 
COMPLEX X(0: L-1) 
REAL A,B 


DO 10 J=K1,K2 

A=REAL(X(J)) 

B=AIMAG(X(J)) 

IF ((A.GE.0). AND. (B.GE.0)) THEN 
R(J-K1)=0 
GO TO 10 

ENDIF 

IF (CA. LE. 0). AND. (B.GE.0)) THEN 


58 





R(J-K1)=1 
GO TO 10 
ENDIF 
IF ((A. LE. 0). AND. ¢B. LE.0)) THEN 
R( JeK1)=2 
GO TO 10 
ENDIF 
R( J-K1)=3 
10 CONTINUE 
RETURN 
END 





Creteterevedevedeve ede de deter re Fete Bee Fede Ne TET NE TEE TENE TENE TE ETE ICSE TENE TEE FETC ETE TCT derksedededededededededetedevede ak 


SUBROUTINE WHAT32(X,N,L,QAMRE, QAMIM,R, PTR, MAX) 


C SUBROUTINE TO FIND CLOSEST QAM CONSTELLATION ELEMENT BY SUBTRACTING 
Cc THE MAGNITUDE OF THE RECEIVED COMPLEX SIGNAL FROM THE PRODUCT OF 
C THE COMPLEX CONJUGATE OF THE RECEIVED SIGNAL AND EACH ELEMENT OF 
C THE QAM CONSTELLATION, AND TAKING THAT CONSTELLATION ELEMENT 

Cc WITH THE SMALLEST RESULT AS AN ESTIMATE OF THE RECEIVED SIGNAL. 
Cc (ONLY USED WITH 32-QAM) 

C 


Cc oe nee a eee ee ee a ee eee ee eT oe ee a ee ee ee eT 
INTEGER N,L,R(0: MAX-1),WHO, PTR 
COMPLEX X(0: L-1) 
REAL NEW,A,NEWR,NEWI, QAMRE( 0: N) , QAMIM(0: N) 
DO 10 I=PTR,PTR+MAX-1 
NEWR=ABS( QAMRE( 0) -REALC(X(I)) 72 
NEWI=ABS( QAMIM(0) -AIMAG(X( 1) ) #2 
NEW=SQRT(NEWR+NEW I ) 
WHO=0 
DO 20 J=1,N 
A=ABS(QAMRE( J) -REAL(X(I)) )#*2 
A=A + ABS(QAMIM( J) -AIMAGCX(I)) )**2 
IF (SQRT(A).LT. NEW) THEN 
NEW= A 
WHO=J 
ENDIF 
20 CONTINUE 


RC I-PTR)=wh 
10 CONTINUE 
RETURN 
END 


Creme vedere sede severe dere tote rete ie vege Fe Be Be He FOIE Te TE TE TET I Te ENE TE IE TE YE VOTE NE TENE TE FEE HED VE BEDE TENE TEE DEVE VE VE DE VE NE TENE TENE VE LENE LE 


SUBROUTINE SYNDRO(R,S,N,NEAREN, OOPS ,GFA, ZEROS , MAX , CODE ) 


C THIS SUBROUTINE COMPUTES THE SYNDROME VECTOR (FOR DECODING) 

C 

C Ditilad ahetannaes SYNDROME VECTOR 

C ZEROS...... NUM™ 2 OF SYNDROME ELEMENTS EQUAL TO ZERO 

CG NEAREN..... FLAG FROM DECODING SUBROUTINE, =1 IF THIS IS THE 
C SECOND SYNDROME COMPUTATION FOR A GIVEN CODEWORD 


Ay 
Cc 





(RECHECKING CORRECTIONS MADE FOR HIGHER ORDER ERRORS) 
OOPS....... FLAG TO DECODING SYNDROME THAT AN ELEMENT OF SYNDROME 
IS NOT EQUAL TO ZERO (WHILE PERFORMING SECOND SYNDROME 
COMPUTATION) 
LOW. G6. ee FIRST ELEMENT OF SYNDROME (DIFFERS BETWEEN LENGTH 31 
AND LENGTH 15 CODES) . 


AAAAAAANA 


INTEGER N, ZEROS ,NEAREN, OOPS ,X,ALPHA, TEMP, SUM, CODE : 
INTEGER R(0: MAX-1),S(-2:2),GFA(0:N,0:N),LOW 
OOPS=0 
DO 5 I=-2,2 
$(1)=0 
CONTINUE 
ZEROS=0 
IF (CODE.GT.2) THEN 
LOW=0 
ELSE 
LOw=-2 
ENDIF 
DO 20 ALPHA=LOW, 2 
SUM=0 
DO 30 I=0,MAX-1 
IF (R(I). EQ. 0) THEN 
GO TO 30 
ENDIF 
TEMP=0 
IF (ALPHA. LT. 0) THEN 
TEMP=ALPHA+N 
X=R(1)-1+TEMP*1 
ELSE 
X=R(1)-1+ALPHA*I : 
ENDIF 
X=MOD(X,N) 
X=X+1 
SUM=GFA(SUM, X) 
30 CONTINUE 
S( Alc HA)=SUM 
IF ((SUM.NE. 0). AND. (NEAREN. EQ. 1)) THEN 
OOPS=1 


ur 


IF (SCALPHA).EQ.0) THEN 
ZEROS=ZEROS+1 
ENDIF 
20 CONTINUE 
21 RETURN 
END 


SUBROUTINE DECTWO(R,S,V,DHAT,E, ERRORS , UNOERR,DUOERR, TRIERR 
+  \FLXED,B,C,K,ERASE,TRACE,N,GFA,GAMMA1, GAMMA2,GAMMA3 
+  ,RUB,NLESSK, ZEROS , CODE , MAX) : 


C DECODING SUBROUTINE FOR DEC/TED CODES 





MAANAAANANAAANAANAAIANANANA 





Reston arid RECEIVED CODEWORD VECTOR 

Dstt chee Cady SYNDROME VECTOR 

Ericecee abit e sy dcane ESTIMATED ERROR VECTOR 

VV 9 Muss to aecieg uhh ESTIMATED CODEWORD VECTOR 

DHAT......... ESTIMATED DATA VECTOR (FROM V) 

ERRORS....... PLAG FOR PARAMETERS OUTPUT IN CASE THERE ARE 
JECODING ERRORS (ERRORS=1) 

UNOERR....... SINGLE ERRORS 

DUOERR....... DOUBLE ERRORS 

TRIERR....... UNCORRECTABLE TRIPLE ERRORS 

ETRED® «aise 840% CORRECTED TRIPLE BURST ERRORS 

Bog Ke,& o.ere% ee PARAMETERS USED IN DOUBLE ERROR CORRECTION 

GAMMA1-3..... PARAMETERS USED IN DOUBLE ERROR CORRECTION 

TRACE. ....... TRACE OF AN ELEMENT OF GF(2***M) 

ERASE........ NUMBER OF ERASURES (IF ERASING TRIPLE-ERRORS IS 


CHOSEN (IF RUB=1) 


wow ewe mete ewww wwe ew ewe mw we ww ee em ww Om weet ewe www Ome BB eM BM BDO wwe HMB wwe ee ee ee 


INTEGER GAMMA1,GAMMA2 ,GANMA3 , UNOERR , DUOERR, TRIERR ,N,MAX 
INTEGER FIXED,NOP,B,C,K,W1,W2,W3,X1,¥2,X3,Y1,Y2,Y3, TEST 
INTEGER ROOT(2),TRACE(0:N),GFA(0:N,0:%!) 

INTEGER R(0:N-1),S(-2:2),V(0: N-1) ,DEAI( 0: N-1) ,E(0: N-1) 
INTEGER Z23,W,X,¥,ERASE,A1,A2, ERRORS, OOPS, ZEROS , RUB, CODE 


eee eee ee eee ee ee eee 


CHECK FOR 3 OR MORE ERRORS FROM SYNDROME... 
IF (ZEROS. GT. 2) THEN 
TRIERR=TRIERR+1 
CALL RUBOUT(RUB ,R,V,DHAT,NLESSK, ZEROS , MAX) 
GO TO 333 
ENDIF 


36 NOP=0 
IF (S(1)*S(-2).EQ.0) THEN 
A1=0 
ELSE 


IF (S(-1)*S(0).EQ.0) THEN 
A2=0 
ELSE 
A2=MOD(S(-1)+8(0)-2,N)+1 

ENDIF 

GAMMA1=GFA(A1,A2) 


IF ((S(2)*S(-2)).EQ.0) THEN 
A1=0 
ELSE 
A1=MOD(S(2)+S(-2)-2,N)+1 
ENDIF 
IF (S(0).E%.0) THEN 
A2=0 * 
ELSE 


A2=MOD(S(0)*2-2,N)+1 
ENDIF 
GAMMA2=GFA(A1,A2) 


IF ((S(0)*S(1)).EQ.0) THEN 
A1=0 


ELSE 
A1=MOD(S(0)+S(1)-2,N)+1 

ENDIF 

IF ((S(2)*S(-1)).EQ. 0) THEN 
A2=0 

ELSE 
A2=MOD(S(2)+S(-1)-2,N)+1 

ENDIF 

GAMMA3=GFA(A1,A2) 


C STEP#4: CHECK FOR A SINGLE ERROR... 


IF ((GAMMA1. EO. 0). AND. (GAMMA2. EQ. 0). AND. (GAMMA3. EQ. 0)) THEN 
Y=MOD(S(1)+N-S(0) ,N) 
E(Y)=S(0) 
UNGERR=UNOERR+1 
GO TO 59 
ENDIF 


C STEP#5: CHECK FOR 3 OR MORE ERRORS... 


a (CCGATIA1. NE. 0). OR. (GAMMA2. NE. 0). OR. (GAMMA3. NE. 0) ). AND. (GAMMA1 
*GAMMA2Q*GANMA3. EQ. 0)) THEN 
TRIERR=TRIERR+1 
ERASE=ERASE+1 
CALL RUEOUT(RUB,R,V,DHAT,NLESSK, ZEROS , MAX) 
GO TO 333 
ENDIF 


C STEP#6: CALCULATE CONSTANTS IN QUADRATIC FORMULA. . 


IF (GAMMA2. EQ. 0) THEN 
B=0 
ELSE 
B=MOD(GAMMA2+N-GAMMA1,N)+1 
ENDIF 
IF (GAMMA3. EQ. 0) THEN 
C=0 
ELSE 
C=MOD(GAMMA3+N-GAMMA1,N)+1 
ENDIF 
IF (C.EQ.0) THEN 
K=0 
GO TO 37 
ENDIF 
IF (B.EQ.1) THEN 


“= 
GO TO 37 
ENDIF 
.=MOD( C-1+3*N-2*(B-1),N)+1 
37 NOP=0 
IF ((TRACE(K). EQ. 1). AND. (K. NE. 6)) THEN 
TRIERR=TRIERR+1 { 
ERASE=ERASE+1 
CALL RUBOUT(RUB,R,V,DHAT,NLESSK, ZEROS , MAX) 
GO TO 333 
ENDIF 


€ LOOK FOR A CORRECTABLE TRIPLE BURST ERROR... 


IF (K.EQ.6) THEN 
Wi=S(-1)-1 
W2=S(1)-1 
/3=MOD(W1+W2,N) 
X1=8(-2)-1 
X2=8(2)+1 
X3=MOD(X1+X2,N) 
Y1=S(0)-1 
IF (N.EQ.15) THEN 
Y¥2=MOD(2**Y143,N) 
Y3=MOD(2**Y1+6 ,N) 
ELSE 
Y2=MOD(2**Y14+20,N) 
Y3=MOD( 2*¥1+23,N) 
ENDIF 
IF ((W3. EQ. Y2). AND. (X3. EQ. Y3)) THEN 
FIXED=FIXED+1 
IF (N.EQ. 15) THEN 
Z3=410D(W2-Y1+3%N-10,N) 
ELSE 
Z3=MOD(W2-Y1+3%N-11,N) 
ENDIF 
E(Z3)=S(0) 
E(23+1)=S8(0) 
E(23+2)=S(0) 
GO TO 59 
ENDIF 


C FIND THE ERRORS USING THE QUADRATIC FORMULA... 


TEST=0 
ROOT( 1)=0 
ROOT( 2)=0 
DO 50 I=1,N 
W=MOD(2%*(I-1),N)+1 
IF (B.EQ.1) THEN 
X=I 
GO TO 41 


ENDIF 
IF (B.EQ.0) THEN 
X=0 
GO TO 41 e 
ENDIF 
X=MOD(B+I-24+N,N)+1 
41 Y=GFA(W, X) 
Z=GFACY,C) 
IF ((2Z.EQ. 0). AND. (TEST. EQ. 0)) THEN 
ROOT(1)=I-1 
42 TEST=1 
GO TO 50 
ENDIF 
IF ((Z. EQ. 0). AND. (TEST. EQ. 1)) THEN 
ROOT(2)=I-1 
43 GO TO 55 
ENDIF 
50 CONTINUE 


55 NOP=0 


C ROOT(1) AND ROOT(2) GIVE DOUBLE ERROR LOCATIONS, NOW FIND MAGNITUDE 
IF (S(0).EQ.0) THEN 
A1=0 
ELSE 
A1=MOD(S(0)-1+ROOT(2),N)+1 
ENDIF 
IF (S(1).EQ.0) THEN 
A2=Al 
ELSE 
A2=GFA(A1,S(1)) 
ENDIF ‘ 
X=MOD(A2+N-B,N)+1 


E(ROOT( 1) )=x 

E(ROOT( 2) )=GFA(S(0) ,X) 

IF (E(ROOT(1))**E(ROOT(2)).EQ. 0) THEN. 
TRIERR=TRIERR+1 
ERASE=EhsSE+1 
CALL RUBOUT(RUB,R,V,DHAT,NLESSK, ZEROS , MAX) 
GO TO 333 

ENDIF 

DUOERR=DUOERR+1 

59 NOP=0 


C STEP#7: NOW ADD ESTIMATED ERROR VECTOR TO RECEIVED VECTOR... 


ERRORS=0 
72 DO 60 I=0,N-1 
V(I)=GFA(R(I) ,E(T)) 
IF (1.GT.4) THEN 
DHAT( 1-5)=V(1) 
ENDIF 
60 CONTINUE 








C STEP##8: MAKE LAST CHECK FOR 3 OR MORE ERRORS... 


NEAREN=1 
ZEROS=0 
OOPS=C6 
CALL SYNDRO(V,S,N,NEAREN, OOPS ,GFA, ZEROS , MAX , CODE) 
IF (OOPS. EQ. 1) THEN 
TRIERR=TRIERR+1 
ERASE=ERASE+1 
DUOERR=DUOERR-1 
CALL RUBOUT( RUB ,R,V,DHAT,NLESSK, ZEROS , MAX) 


GO TO 333 
ENDIF 
333. RETURN 

ane 

Cisdedevedsvede dey: ‘sdevetededededevedevedess setedededkedeveve severe te sewed Seledesededededetederevededetetsielesedeksevetedesetede 
SUBROUTINE OECORE CS: V, ZEROS, RUB, R, pie N,NLESSK,MAX,GFA, CODE, 

+ UNOERR DUCERR) 
C 
C DECODING SUBROUTINE FOR SEC/DED CODES 
Cc 
C Riad aest tied RECEIVED CODEWORD VECTOR 
C Soar etebeud aeheiardes SYNDROME VECTOR 
C Vi ahand soe Ses, 808 ESTIMATED CODEWORD VECTOR 
C DHAT......... ESTIMATED DATA VECTOR (FROM VY) 
C DECODING ERRORS (ERRORS=1) 
C UNDERR. .. 2... SINGLE ERRORS 
C DUOERR....... DOUBLE ERRORS 
C CHOSEN (IF RUB=1) 
Cc 
Cc same eeeeeas ese e wen mmmew ere newer seer ewe eceseat ee wm ewes s wre wrenawreweeweeevawnecacae 


INTEGER ZEROS ,N,DUOERR,X,Y,UNOERR,MAX,CODE,RUB,NLESSK 
INTEGER S(-2 ),GFA(0O:N,0:N),V(0: MAX-1) ,DHAT(0: MAX-1) ,R(O: MAX-1) 
INTEGER NEAREN, OOPS 


IF (ZEROS. EQ. 1) THEN 
DUOERR=DUOERR+1 
CALL RUBOUT(RUB,R,V,DHAT,NLESSK, ZEROS , MAX) 
GO TO 150 
ENDIF 
S(-2)=0 
S(-1)=0 
DO 10 I=0,MAX-1 
V(I)=R(I) 
10 CONTINUE 
IF ((S(0).NE. 0). AND. (ZEROS. EQ. 2)) THEN 
V(0)=GFA(S(0) ,RC0) ) 
UNOERR=UNOERR+1 
GO TO 100 
ENDIF 
IF ((S(1).NE. 0). AND. (ZEROS. EQ. 2)) THEN 
VC 1)=GFA(S(1),RC1)) 
UNOERR=UNOERR+1 


GO TO 100 
ENDIF 
IF ((S(2).NE. 0). AND. (ZEROS. EQ. 2)) THEN 
V(2)=GFA(S(2) ,RC2)) 
UNOERR=UNOERR#+1 
GO TO 100 
ENDIF ‘ 
X=MOD(S(1)+N-S(0),N) 
Y=MOD(S(2)+N-S(1),N) 
IF (X.EQ.Y) THEN 
C X=X+3 
V(X)=GFA(R(X) ,8(0) ) 
UNOERR=UNOERR+1 
GO TO 100 
ELSE 
DUOERR=DUOERR+1 
CALL RUBOUT( RUB,R,V,DHAT,NLESSK, ZEROS , MAX) 
GO TO 150 
ENDIF 
C LAST CHANCE TO CHECK FOR 2 OR MORE ERRORS... 
100 NEAREN=1 
ZEROS=0 
OOPS=0 
CALL SYNDRO(V,S,N,NEAREN , OOPS ,GFA, ZEROS , MAX , CODE) 
IF (OOPS.EQ.1) THEN 
DUOERR=DUOERR+1 
UNOERR=UNOERR-1 
CALL RUBOUT(RUB,R,V,DHAT,NLESSK, ZEROS , MAX) 
GO TO 150 
ENDIF 
DO 20 I=3,MAX-1 
DHAT( I-3)=V(1) 
20 CONTINUE 
150 RETURN : 
END 


Cretevededededs totes: devedevededededecedevevedededete dete de ete Te see Te Te See Ne TEER TESE TLE DOSE SESE SE SEVEN TENT E TOD ELE DODD NCTE TERRE 


SUBROUTINE RUBOUT(RUB,R,V, DHAT, NLESSK, ZEROS , MAX) 


C THIS SUBROUTINE ERASES CODEWORDS WITH MORE ERRORS THAN A GIVEN CODE 
C CAN CORRECT. 
C 
Cow ccc ce ence cc ecw wenn een nnn eee nnn ee een nnn came enna ee coese 
INTEGER NLESSK, ZEROS , MAX 
INTEGER RUB,R(O: MAX-1) ,V(0: MAX~1) , DHAT(0: MAX-1) 
IF ((RUB. EQ. 0). OR. (ZEROS. EQ. NLESSK)) THEN 
DO 14 I=0,MAX-1 
VC I)=R(T) 
IF (1I.GT.NLESSK-1) THEN ji 
DHAT( I-NLESSK)=R(1) 
ENDIF 
14 CONTINUE , 


GO TO 50 





ENDIF 
IF (RUB. EQ. 1) THEN 
DO 101 I=0,MAX-1 
V(I)=-1 
IF (I. LE. MAX-NLESSK-1) THEN 
DHAT(I)=-1 
ENDIF 
101 CONTINUE 
ENDIF 
50 RETURN 
END 


ferloovlevlevlovledeniedovevlovevlevlevovedevleslodesledeslerlevededadededededeviovloleyevevesesloviodedesovededesesiedeveves. ledlovedededevedeseveslevedesiedeve 
Veee ese Tere sete Pete te eee LTTE TELL EE TEEN VEDE TEESE DEEL LELELE ELEN SESE ELE EN DEEPA AED AEBS AEE 


SUBROUTINE OUTPUT(S ,GAMMA1 ,GAMMA2 ,GAMMA3,R,V,B,C,K,CW,D, 
+ BAUDNO , DHAT, BDTYPE , MAX) 


FOR BYTE ERRORS, WRITE PERTINENT DATA TO RSOUTPUT FILE... 


weew ewe emr eww ew men ee eee mmm eee eee we meme e wm wwe ewe we wn we ewe meme ne weer e eee ema ene 


INTEGER GAMMA1,GAMMA2 ,GAMMA3,B,C,K,MAX 
INTEGER S(-2:2),R(0: MAX) ,V(0: MAX) , CWO: MAX) 
INTEGER DCO: MAX) , BAUDNO, DHAT( 0: MAX) , BDTYPE 


WRITE( 6%)! oon rena n nen cence een nen eee eee een ene en ee 


WRITE(6,*)'BAUDTYPE=' , BDTYPE 
WRITE(6,*)} ‘CODEWORD NUMBER=' , BAUDNO 
WRITE(6 , 904) 
904 FORMAT( 'SYNDROMES' ) 
WRITE(6,905) S(-2),S(-1),S(0),S(1),8(2) 
905 FORMAT(1X,12,3X,12,3X,12,3X,12,3X,12) 
WRITE(6,*)' 
WRITE(6,*)' ' 
WRITE(6,906) 


906 FORMAT( GAMMAS' ) 
WRITE(6,907) GANMA1,GAMMA2 ,GAMMA3 
907 FORMATC SS 12 5X 12,524 22) 


WRITE(6,*)'B,C,K' 
WRITE(6,*)B,C,K 
WRITE(6,*)' ' 
WRITE(6,908) 
908 FORMAT(11X,‘'I',10X,'DATA',7X,'ENCODE',7X,'RECD',8X,'ESTCW', 
+ 5X, 'ESTDATA') 
DO 71 J=0,MAX-1 
WRITE (6,*) J,D(J),CW(J),R(J),V(J), 
+ : DHAT( J) 
71 CONTINUE 
RETURN 
END 


Crevederevedetede decade dec dete rece tote Bete LENT ee TOL eL 8 BENE LELE LE SE LET ELE SESE SELEVESE PELE LE LESSEN DEED SER IVOLEE LINN 
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SUBROUTINE INTRAN(D1,D,MAX,NLESSK) 


SUBROUTINE TO TRANSLATE REAL DATA INTO INTEGERS... 
(SINCE CODER/DECODER DEALS ONLY WITH INTEGERS) 


Br ansceey eaten INTEGER OUTPUT 
Dlsc.ch.os aes REAL INPUT 


AAQARAAANaA 
oO 
oO 


wens eeen nme nn eee ee see ewer eww mn ene ew ee Feo ee wwe tee we eeon se emm ete ew eee wee e se = 


INTEGER MAX,NLESSK 
INTEGER D(0: MAX-NLESSK-1) 
REAL D1(0: MAX-NLESSK-1) 


DO 1 I=0,MAX-NLESSK-1 
DO 2 J=0,MAX-1 
IF (D1(1I).EQ.(J/1.0)} THEN 


DC 1)=J 
ENDIF 
2 CONTINUE 
1 CONTINUE 
RETURN 
END 


Cistevededetedeve decade vets fete te te Te Fe TeV Te eee Te TENET BE ETE TELE TE TELE TE LEVEE VE TE TEDE NE VED NE TELE DE TELE LE TELE NE NNDB LCN ICT NE 


SUBROUTINE BITREV(N,M,XTMP,X) 


C THIS SUBROUTINE ACCEPTS AS INPUT THE COMPLEX ARRAY 
C CONTAINING THE VALUES XTMP() WHICH WERE READ FROM 
C THE INPUT FILE. THE OUTPUT OF THIS SUBROUTINE IS 
C THE COMPLEX ARRAY X() WHICH CONTAINS THE INPUT 
C VALUES IN 'BIT-REVERSED' ORDER. 
C (FROM THE LIBRARY OF PROFESSOR STiUM AT NPS) 
Cc 
CT TCL 
COMPLEX XTMP(0:N-1),X(0:N-1) 
DO 10 K=0,N-1 
NEWADR=0 
MADDR=K 
DO 20 I=0,M-1 
LRMNDR=MOD( MADDR, 2) 
NEWADR=NEWADR+LRMNDR* 2****(M-1-1) 
MADDR=MADDR/2 
20 CONTINUE 


XC(NEWADR )=XTMP(K) 
10 CONTINUE 


RETURN 
END 
’ 
Crevederedeveve levered terete re tec ere Le Lee Te eR ORL ETE EEE LER EEL EPO VCTE LE TELEVE DER ELEN EN NEES SER SII 


SUBROUTINE FFT(N,M,X) 
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THIS SUBROUTINE ACCEPTS AS INPUT THE COMPLEX ARRAY X(), 
CALCULATES THE FAST FOURIER TRANSFORM (FFT) OF THE 
ARRAY, AND RETURNS THE RESULTING SEQUENCE IN THE 
ORIGINAL ARRAY X(). 

(FROM THE LIBRARY OF PROFESSOR STRUM AT NPS) 


AAgAAaAaAN 


wseaemawereeremrzweee een wesw este en eww ewe eee ewe ww ww mw we ww we mw wom et wwe we ote eee te 


COMPLEX X(0:N-1),W, TMP 
PI=4. O*ATAN( 1. 0) 
EN=N 


DO 50 L=1,M 
ISPACE=2***L, 
S=N/ISPACE 
IWIDTH=ISPACE/2 
DO 40 J=0,C(IWIDTH-1) 
R=S*J 
ALPHA=2. 0*PI*R/EN 
W=CMPLX(COS( ALPHA) , -SIN(C ALPHA) ) 
DO 30 ITOP=J,N-2,ISPACE 
IBOT=ITOP+IWIDTH 
TMP=XCIBOT)*W 
XC IBOT)=XC ITOP)-TMP 
XC ITOP)=KCITOP)+TMP 
30 CONTINUE ; 
40 CONTINUE 
50 CONTINUE 
RETURN 
END 
Crsveteveds rede deve de ted ete ce rete Tobe tote Te VEG ETE TE TE TENE TE TE TE TE TENE TELE TEE VE TELE TE TE TE VE TE NEVE SET EVEN ELON LEE LESTE TEE LE ETI NE PONE 
SUBRGUTINE INVFFTCN,M,X) 


THIS SUBROUTINE ACCEPTS AS INPUT THE COMPLEX ARRAY 
X() AND RETURNS THE INVERSE FFT OF THE ARRAY. 
(FROM THE LIBRARY OF PROFESSOR STRUM AT NPS) 


aANAN|ANANA 


COMPLEX X(0:N-1) 
EN=N 


C CALCULATE THE COMPLEX CONJUGATE OF THE INPUT DATA. 
DO 70 I=0,N-1 
XC I)=CONIGCX(1)) 
70 CONTINUE 
C CALCULATE THE FAST FOURIER TRANSFORM OF THE ARRAY. 
CALL FFT(N,M,X) 


C CALCULATE THE COMPLEX CONJUGATE OF THE RESULTING ARRAY. 
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DO 80 I=0,N-1 
XCI)=CONJG(X(1))/EN 
80 CONTINUE 


RETURN 
END { 
Ciedededededede detects tee Tete te Fe te Fe Fe TNE FETE LENE ETE TENE TEE TE TE ETE ETE TENE TERETE TECTED TE NE DETTE VETERE SELENE TERI RICK 


SUBROUTINE GAUSS( AVG, VAR, PTS ,AWGN, TEMP1,TEMP2) 


C€  GAUSS-SUBROUTINE TO OUTPUT A GAUSSIAN RANDOM VARIABLE ARRAY WITH 
Cc MEAN=AVG AND VARIANCE=VAR 

C 

Cc AVG....... MEAN OF THE GAUSSIAN RANDOM VARIABLE 

Cc VAR. ...... VARIANCE OF THE GAUSSIAN RANDOM VARIABLE 

C 

C 


INTEGER PTS 

REAL AVG,VAR,TEMP1(0: PTS-1) , TEMP2(0: PTS-1) 
REAL A,B,PI 

REAL AwWGN(0: PTS-1) 


PI=4. O*ATAN(1. 0) 
CALL RNUN(PTS,TEMP1) 
CALL RNUN( PTS, TEMP2) a 
DO 10 I=0,PTS-1 
A=-2. 0°**ALOG(TEMP1(I)) 
B=COS( 2**PI*TEMP2(1)) 
AWGN( I)=SQRT(A)**B**SQRT( VAR)+AVG 
10 CONTINUE 
RETURN 
END 





evevevededededevevede ve ve de ve dete ce tesa te Te Ree MeN ese Ne Te eK Te TeV He Tee eM TE eM Sete se Ne Ne THe 


evededevedededesedeve dete ve deve deve de sess Ne MeV eT Te Ne TTT LEN TTT TE 
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